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Low-frequency  (LF)  seabed  scattering  at  low  grazing  angles  (LGA)  is  almost  impossible  to  directly 
measure  in  shallow  water  (SW),  except  through  inversion  from  reverberation.  The  energy  flux 
method  for  SW  reverberation  is  briefly  introduced  in  this  paper.  The  closed-form  expressions  of 
reverberation  in  an  isovelocity  waveguide,  derived  from  this  method,  indicate  that  in  the  three- 
halves  law  range  interval  multimode/ray  sea  bottom  scattering  with  different  incident  and  scattering 
angles  in  forming  the  reverberation  may  equivalently  be  represented  by  the  bottom  backscattering 
at  a  single  range-dependent  angle.  This  equivalent  relationship  is  used  to  derive  the  bottom  back- 
scattering  strength  (BBS)  as  a  function  of  angle  and  frequency.  The  LF&LGA  BBS  is  derived  in  a 
frequency  band  of  200-2500  Hz  and  in  a  grazing  angle  range  of  1.1°-14.0°  from  reverberation 
measurements  at  three  sites  with  sandy  bottoms.  This  is  based  on  three  previous  works:  (1)  The 
closed-form  expressions  of  SW  reverberation  [Zhou,  (Chinese)  Acta  Acustica  5,  86-99  (1980)];  (2) 
the  effective  geo-acoustic  model  of  sandy  bottoms  that  follows  the  Biot  model  [Zhou  et  al.,  J. 
Acoust.  Soc.  Am.  125,  2847-2866  (2009)]  and  (3)  A  quality  database  of  wideband  reverberation 
level  normalized  to  source  level  [Zhou  and  Zhang,  IEEE  J.  Oceanic  Eng.  30,  832-842  (2005)]. 
©2012  Acoustical  Society  of  America.  [http://dx.doi.org/10.1121/L3693645] 
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I.  INTRODUCTION 

Acoustic  reverberation  in  shallow  water  (SW)  wave¬ 
guides  has  become  a  hot  research  topic  in  recent  years.  It 
involves  a  complex  process  of  both  two-way  sound  propaga¬ 
tions  and  bottom  scattering.  Figure  1  shows  a  flow  chart  of 
seabed-dominated  SW  reverberation.  The  seabottom  geoa¬ 
coustic  parameters  and  velocity  profiles  in  the  water  column 
control  two-way  propagations  that  can  now  be  dealt  with  by 
using  existing  theoretical  formulae  and  numerical  methods, 
as  summarized  in  a  monograph  by  Jensen  et  al. 1  Thus,  as  the 
main  source  of  SW  reverberation,  seabottom  scattering 
becomes  a  kernel  problem  in  the  development  of  reverbera¬ 
tion  models. 

In  1999,  Preston  indicated2  “Much  progress  has  been 
made  in  our  understanding  of  reverberation.  However,  there 
remain  important  unanswered  questions  and  a  real  scarcity 
of  high-quality  basic  research  data  sets.”  Figure  2  of  Ref.  2 
shows  that  currently  available  BBS  data  at  LF&LGA  vary  as 
much  as  30  dB  at  a  given  angle.  This  large  variation  might 
be  caused  by  different  sediment  properties,  but  it  might  also 
be  due  to  questionable  reverberation  data,  due  to  the  theoret¬ 
ical  formula  of  reverberation,  or  the  seabed  geoacoustic 
models  used  in  some  scattering  inversions. 

The  derivation  of  a  reliable  long-range  LF  reverberation 
model  requires  a  reliable  LF&LGA  bottom  scattering  model 
that  is  proven  by  experimental  data.  The  LF&LGA  seabed 
scattering  is  almost  impossible  to  directly  measure  in  shallow 
water,  and  is  generally  derived  from  mid-  and  long-range 
reverberation.  As  Eller  argued,  “Unless  the  assumptions  used 
in  extracting  the  scattering  coefficients  from  the  received 
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reverberation  are  very  carefully  chosen,  the  resulting  coeffi¬ 
cients  will  not  represent  valid  environmental  parameters,  in¬ 
dependent  of  how  the  data  were  obtained  and  analyzed.”3 

In  the  past  30  years,  one  of  the  major  accomplishments 
in  ocean  acoustics  has  been  our  understanding  of  seabed 
scattering.  This  has  been  summarized  in  a  monograph  of  the 
ONR  Underwater  Acoustics  Series,  entitled  “ High- 
Frequency  Seafloor  Acoustics.”4  As  its  introduction  indi¬ 
cates,  this  monograph  emphasizes  “high  frequencies,  very 
roughly,  frequencies  from  10  kHz  to  1  MHz.”  A  natural 
question  is  raised  from  this  book:  Are  those  seabed  scatter¬ 
ing  models,  which  are  well  developed  for  seabed  boundary 
roughness  and  sediment  inhomogeneity,  suitable/applicable 
to  low  grazing  angles  and  low  frequencies  of  100-3000  Hz? 
To  answer  this  question,  we  also  need  the  LF  &LGA  seabed 
scattering  data  that  are  currently  in  “a  real  scarcity.”2 

In  addition  to  a  reliable  reverberation  model,  there  are 
two  other  absolute  requirements  for  the  inversion  of  the 
LF&LGA  scattering  from  reverberation:  (1)  carefully  cali¬ 
brated  reverberation  measurements  in  a  broad  band  and  (2)  a 
ground  truth  about  the  seabed  geo-acoustic  model  in  the 
same  sea  area.  Poorly  calibrated  reverberation  data  result  in 
incorrect  values  of  the  seabed  scattering.  A  previous  paper 
showed  that  a  decrease  (increase)  in  bottom  reflection  loss 
can  be  compensated  for  by  increasing  (decreasing)  the  angu¬ 
lar  index  of  bottom  scattering.5  That  is,  there  is  an  uncer¬ 
tainty  caused  by  an  internal  coupling  between  a  seabed 
geoacoustic  model  and  a  scattering  model. 

Taking  the  above-mentioned  three  essential  conditions 
of  reverberation  inversions  into  account,  this  paper  derives 
the  LF&LGA  seabed  scattering  strength  vs  angle  and 
frequency  in  SW  with  sandy  bottoms,  based  on  three  previ¬ 
ous  works:  (1)  The  analytical  (closed-form)  expressions  of 
SW  reverberation;6-8  (2)  an  effective  geo-acoustic  model 
of  sandy  bottoms  that  follows  the  Biot  model  as  well  as 
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FIG.  1 .  (Color  online)  Flow  chart  of 
seabed-dominated  reverberation  in 
shallow  water. 


other  physics  based  models;9’10  and  (3)  a  quality  database 
of  wideband  reverberation  level  (RL)  normalized  to  source 
level  (SL).11  This  paper  is  organized  as  follows:  Sec.  II 
briefly  introduces  the  angular  spectrum  (energy  flux) 
method  for  SW  reverberation  and  the  closed-form  expres¬ 
sions  of  reverberation  in  isovelocity  shallow  water.  A  theo¬ 
retical  model  and  a  methodology  which  is  used  to  derive 
the  LF&LGA  bottom  backscattering  strength  are  described 
in  detail  in  Sec.  III.  Section  IV  describes  the  seabed  geoa¬ 
coustic  model  used  in  the  derivation  of  BBS.  The 
reverberation-derived  LF&LGA  BBSs  as  a  function  of  fre¬ 
quency  and  angles  for  three  sites  with  sandy  bottoms  are 
given  in  Sec.  V.  Finally,  Sec.  VI  summarizes  and  discusses 
the  results  of  this  paper. 


II.  THE  CLOSED-FORM  EXPRESSIONS  FOR  OCEAN 
REVERBERATION  INTENSITY  IN  SHALLOW  WATER 

The  angular  spectrum  (energy  flux)  method  for  model¬ 
ing  SW  reverberation,  based  on  the  WKB  approximation  of 
normal-mode  theory,  was  first  presented  three  decades  ago 
in  archived  Chinese  journals.6,7’12'13  These  included  closed 
form  expressions  for  reverberation  in  SW  with  a  down 
refracting  and  isovelocity  profile  or  with  a  varied- 
depth. ^ 6_8- 13-16  £)ue  to  its  simplicity  and  intuitive  nature,  the 
energy  flux  method  has  further  been  developed  and  widely 
used  for  SW  reverberation  calculations, 17-30  mainly  due  to 
works  by  Harrison,  Holland  and  Ainslie. 18,20  21  23-30  All  the 
closed-form  expressions  in  isovelocity  shallow  water, 
derived  by  Zhou,6  Harrison,18  Ainslie27,30  and  Holland,24  are 
almost  identical,  having  only  minor  differences  in  constant 
coefficients.  In  this  paper,  Zhou’s  closed-form  expressions 
for  reverberation  in  isovelocity  water  are  used  to  extract  the 
frequency  and  angular  dependence  of  the  BBS  at  LF&LGA 
from  long-range  reverberation  data. 


approximation  to  the  normal-mode  solution  of  the  wave 
equation,  was  in  Chinese.6  Its  English  translation  has  been 
unavailable  to  all  but  a  very  few  western  colleagues.  It  is 
worthwhile  to  briefly  introduce  this  method  here. 

A  source  level  (SL)  of  0  dB  will  be  assumed  in  all  sub¬ 
sequent  expressions.  The  sound  field  intensity  in  shallow 
water,  normalized  to  SL,  can  be  expressed  as  a  sum  of  nor¬ 
mal  modes: 

7(r,  z;  z0)  =  -  V 

r  I  kn 

y-  y-  (lL(zo)$*,(zo)<l>„(z)0*,(z) 

2-^  2^  ,1/2,  *1/2 

n^m  & m 

X  g-[(Pn+Pm)+Kkn-km)]r 

where  is  the  eigenfunction,  k„  is  the  longitudinal  wave 
number,  [in  is  the  modal  attenuation  factor  of  the  nth  mode, 
zo  is  the  source  depth  and  z  is  the  receiver  depth.  The  first 
term  of  Eq.  (1)  is  an  incoherent  summation  of  normal- 
modes.  It  describes  the  average  intensity  of  the  sound  field. 
Because  SW  environments  are  complex  and  variable  in  both 
space  and  time,  it  is  difficult  to  predict  or  measure  the  fine 
structure  of  a  SW  sound  field  associated  with  modal  interfer¬ 
ence.  Thus,  for  many  practical  applications  such  as  sonar 
performance  predictions,  the  second  term  of  Eq.  (1),  which 
describes  the  interference  fine  structure,  can  be  neglected. 

According  to  the  WKB  approximation,  the  smoothed 
(slowly  varying)  energy  depth  distribution  of  the  nth  mode 
can  approximately  be  expressed  as 


l$„(z)|2 


2 

S„  tan  0„(z)  ’ 


(2) 


Here,  S„  is  the  cycle  distance  of  the  mil  mode-ray  given  by 


A.  General  expressions  of  the  energy  flux  method  for 
average  SW  sound  intensity  and  reverberation 


S„ 


2 


',,n  dz 
tan  9n  (z) 


(3) 


The  earliest  detailed  derivation  of  the  energy  flux 
method  for  SW  reverberation,  derived  from  the  WKB 


and  6n  is  the  mode  grazing  angle.  The  eigenvalue  kn  of  the 
nth  mode  satisfies  the  relationship 
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Jk2{z)  -  k2dz  +  E(n  +  s,hi  =  2/771, 


we  can  derive  a  general  expression  for  the  average  field  in¬ 
tensity  in  shallow  water  given  by 


where  i;n  and  rjn  are  the  upper  and  lower  turning  (or  reflec¬ 
tion)  depth  of  the  nth  mode,  e^n  and  are  the  correspond¬ 
ing  phase  shifts,  and  k(z)  =  c o/c(z).  Neglecting  other 
possible  losses  except  bottom  reflection  loss,  the  modal 
attenuation  factor,  /?„,  can  be  expressed  approximately  in 
terms  of  the  bottom  reflection  coefficient  V(9): 


ln|V(0n)| 
Sn(6n )  +  <5„ 


(5) 


where  5„  is  a  beam  displacement  on  the  bottom  reflection.  If 
the  water  depth  H  satisfies  kH  3>  I ,  Sn  <C  Sn;  i.e.,  when  the 
frequency  is  relatively  high,  dn  is  negligible. 

Each  normal-mode  can  be  decomposed  into  a  pair  of 
up-  and  down-going  local  quasi-plane  waves.  We  shall  call 
the  local  plane  wave  an  equivalent  “mode-ray.”  The  grazing 
angle  of  the  mode-ray  and  the  eigenvalue  of  the  nth  mode 
are  related  by 


kn  =  k{z )  cos  6„(z). 


(6) 


Ignoring  the  minor  effects  of  e^n  and  ,  differentiating  kn 
with  respect  to  n  in  Eq.  (4),  and  using  Eqs.  (3),  (6)  and  the 
identity, 


d 

dn 


•n„ 


/(z,  n)dz  = 


f'(z,n)dz  +f(r\n,ri) 

{„ 


d>1n 

dn 


d_L 

dn 


we  obtain  the  known  relationship:1 


(7) 


dkn  2n 

dn  Sn 


(8) 


Using  (6)  and  (8),  we  have 


Snk(zo)sm6n(zo)d0n(zo) 

dn  = - - - .  (9) 
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The  incoherent  summation  over  modes  can  be  changed  to  an 
integration  over  the  mode-ray  angle  9  .  Proceeding  from  (1), 
(2),  (5),  (8),  and  (9),  and  adding  a  term  of  water  absorption, 


7(r,  z;z0)  =  -e 

 2e^r 
r 


2e21n|V(0)|rA„(fl) 


d9„(z0) 


s„  tan  0„(z) 
lapsi.9,  r ■  z,  z0)d9n (z0) 


(10) 


where 


2e~2Pnr  2e2ln\v(-6^rlsdB) 

laps  {9,  T ,  z;  Zo)  =  — —  .  .  =  — —  .  .  (11) 

Sn  tan  9n  (z)  S„  tan  9n  (z) 

Using  different  methods,  both  Brekhovskikh31,32  and  Smith1’ 
obtained  a  similar  expression.  In  (10)  and  (11),  a  is  the  water 
column  absorption  coefficient;  9(z0 )  is  the  grazing  angle  of 
the  mode-rays  at  the  source  depth  and  9(z)  is  the  grazing 
angle  at  the  receiver  depth.  The  quantity  defined  by  Eq.  (11), 
Iaps  is  referred  to  as  the  “angular  power  spectrum”  of  sound 
propagation.  Smith  also  called  it  the  “energy  density  in  direc¬ 
tion  of  arrival.”3’  The  average  sound  intensity  in  shallow 
water  depends  on  water  absorption  and  two  other  parameters: 
the  bottom  reflection  coefficient  V(9)  and  the  cycle  distance 
S„  of  the  mode-ray.  According  to  Guthrie  and  Tindle,  the 
“mode  rays”  (eigenrays)  associated  with  normal  modes 
are  spatially  fuzzy.  The  “mode-rays”  obey  Snell’s  law: 
kn  =  k{z)  cos  9n  (z)  =  k(zo)  cos  9„(zo).  The  interference  of 
many  adjacent  modes  forms  a  physical  ray  with  limited  width 
that  converges  to  a  geometric  ray  at  high  frequency. ’4 

Equation  (10)  is  a  simple  and  intuitive  expression  for 
calculating  averaged  SW  field  characteristics  such  as  sound 
propagation,  reverberation,  noise  and  their  spatial  coher¬ 
ence.  It  is  only  an  angular  weighting  process.  For  calculat¬ 
ing  the  average  reverberation  intensity  in  shallow  water, 
the  two-way  propagation  angular  spectra  and  the  classic 
seabed  scattering  function  are  well  connected  in  the  angular 
domain. 

For  a  frequency  band  this  paper  is  interested  in,  and  at 
low  order  modes-dominated  long-range,  the  water  absorption 
loss  may  be  approximated  as  the  horizontal  range  times  the 
absorption  coefficient  a.  Using  the  single  scatter  approxima¬ 
tion,  the  average  reverberation  intensity  R  can  be  expressed 
in  terms  of  Iaps  by6 


R(r.z::0)  = 


e  A*Mb(0.*)  — l^.r.z.H)  dddf 


r 


rce 


-•/;  2rHn|V'{fl)|(-/.t„(0) 

/■  S„  x  tan  0„(H) 


xA  xMh[0n(H),cpm(H)}  x 


?~*r  2  e2ln  I  Uv)  !/•/■»*,  (pi 
r  Sm  x  tan  tpni(z) 


dO„U)d<Pm(H) 


(12) 
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where  Iaps(9,  r,  II:  zq)  is  the  incident  angular  power  spectrum 
at  a  distance  of  r,  incident  upon  a  seabed  area  A; 
IapS{(p.  r.  z:  II)  is  the  angular  power  spectrum  returned  to  a 
receiver  at  a  depth  of  z  from  the  same  bottom  scattering 
area;  77  is  the  water  depth;  M/,[9n(H),  cpm(H)}  is  the  bottom 
scattering  function;  0„(H)  is  the  incident  grazing  angle  of 
mode  n  at  the  bottom;  and  < pm(H)  is  the  scattering  angle  of 
mode  m. 


B.  The  average  sound  intensity  in  the  Pekeris 
waveguide  (isovelocity  water) 

The  plane-wave  reflection  coefficient  from  the  half-space 
seabed  at  grazing  angle  9  can  be  approximated  by35-38 

,  f  Q°  0  <0<9C 

nlV^-)l  -  {  -  In  |  Vo  |  =  const  9C  <  9  <  n/ 2,  (13) 

where  9c  is  the  critical  angle,  Q  is  a  seabed  reflection 
loss  factor;  the  seabed  reflection  loss  will  be 
— 20 loglo|V(0) |  =  8.686(30  (dB/rad).  Vo  is  the  bottom 
reflection  coefficient  at  large  grazing  angles  close  to  the 
normal  incidence.  In  isovelocity  shallow-water,  the  cycle 
distance  of  the  nth  mode-ray,  S„  =  277/tan(0„).  From  (10) 
and  (13),  the  average  sound  intensity  as  a  function  of  range 
in  the  isovelocity  shallow-water  can  be  divided  into  four 
regions,  with  different  “decay  laws”  expressed  (omitting 
sound  absorption  in  the  water  column)  by6’35-37 

(a)  For  /'o  <  r  <  r\ ,  spherical  spreading  region  ~  r 

,  ,  ,  4  1 

a(0  ~  ^ |  yQ|  r2  ' 


was  defined  first  by  Brekhovskikh  as  the  “three-halves  law” 
region.16  It  is  roughly  hundreds  of  meters  to  several  tens  of 
thousands  of  meters  from  a  sound  source.37 

C.  The  closed-form  expressions  of  reverberation  in 
the  Pekeris  model  (Ref.  6) 

The  descriptions  of  seabed  scattering  by  Urick39  and 
Mackenzie-Lambert40  were  extended  to  more  general  angu¬ 
lar  dependence  as  follows. 6-8  When  a  plane  wave  with  an  in¬ 
tensity  of  7,  (0)  is  incident  upon  unit  bottom  element  dA  at  a 
grazing  angle  9,  the  backward  scattering  intensity  at  one  me¬ 
ter  from  an  elementary  area  in  the  cp  direction  can  be  written 
as 

dls  =  n h{9>)  sin”’  9 sinA  cpdA .  (19) 

Here  fi  is  the  generalized  Lambert  coefficient,  m  and  k  are 
the  angular  indexes  for  incident  and  scattering  waves, 
respectively. 

The  insonified  area  A  is  defined  as  A  m  nrcx.  Making 
the  small  angle  approximation,  neglecting  angles  greater 
than  the  critical  angle,  and  with  the  aid  of  Eqs.  (13),  (19)  and 
the  definition  of  A,  Eq.  (12)  may  be  solved  to  produce  the 
following  closed-form  solutions,  as  described  in  Eqs.  (4.1)  to 
(4.3)  of  Ref.  6. 

(1)  For  the  Lambert  scattering  model  (m  =  k  —  1): 

Rn{r)=^Le-2^{l_e-(  dr/mf.  (20) 

4  Q-r* 

(14)  (2)  For  the  Omnidirectional  model  (m  =  1  ,  k  =  0): 


(b) 

(c) 

(d) 


For  r\  <  r  <  r 2,  cylindrical  spreading  region  ~  r  1 , 

(15) 


,n  29c  1 


For  /'2  <  r  <  /"3,  three-halves  law  region  ~  r  3/2 
c(,)  y  QH  r3/2  ‘ 

For  r  >  >‘3,  single  mode  (1st  mode)  decay  region 


(16) 


Rw(r) 


^7I3/2CT 

4/7‘/2Q3/2r5/2 

x  e-2”'(l  -  c-(Gr/#O^)0 


(21) 


(3)  For  the  Angle-independent  scattering  model  (m  =  k  =  0): 

lQrr  ' 


*ooM  =  ^V2^2 


41 1  (Jr - 


(22) 


b{r) 


2k  1 

idf-7 


exp 


Qn~r 

k2H3 


(17) 


Three  transition  distances  in  the  Pekeris  model  are 
defined  as 


where  x  is  the  pulse  duration  of  the  transmitting  signal, 
and  c  is  the  sound  speed  in  the  water  column.  The  sub¬ 
script  indices  of  Rmk,  m  and  k,  are  the  angular  index  of 
incidence  and  the  scattering  angle,  respectively.  (f>{x)  is 
the  probability  integral 


77 


r  1  = 


(-ln|V*|)0c 


ri  = 


77 

Q^, 


,2  ’ 


t’3  = 


k2H 3 
Qn2 


(18) 


Hx)  = 


ex 


p(-r)dt. 


where  k  is  the  wave  number.  In  general,  active  sonar  predic¬ 
tions  are  most  interested  in  the  mode  stripping  region  C  that 


As  Sec.  IV  of  Ref.  6  indicates,  “although  the  above- 
mentioned  formulae  are  derived  for  n  =  ( m  4-  k)  =  0,  1  and 
2,  our  numerical  calculations  show  that  the  results  are 
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approximately  applicable  to  arbitrary  values  of  the  bottom 
scattering  index  n,  i.e.,  a  combination  of  any  k  and  m,  inte¬ 
ger  or  decimal.  It  only  requires  a  minor  modification  for  the 
constant  coefficient.”  Thus,  the  closed-form  expressions  of 
long-range  reverberation  in  isovelocity  shallow  water  were 
generalized  to  Eq.  (4.15)  and  Eq.  (4.17)  of  Ref.  6  as 
follows: 

( 1 )  When  /'2  >  r  >  r\ 


R(r) 


UTZCldf^  ^  c_2Xvl.  1 


(23) 


(2)  When  /'3  >  r  >  r 2 


R(r) 


jU7t(2-"/2)  CT 
4Q(l+n/2)H(l-n/2)e 


1 

r(2+"/2) ' 


(24) 


Here  the  bottom  scattering  index  n  =  m  +  k,  i.e.,  it  is  a 
combination  of  any  k  and  m,  integer  or  decimal.  Harrison,18 
Ainslie27  and  Holland24  have  derived  expressions  for  an 
isovelocity  SW  that  are  similar  to  Eq.  (23)  and  Eq.  (24) 
with  a  minor  difference  in  the  constant  coefficients. 


Using  Eq.  (16)  and  Eq.  (26),  Eq.  (24)  may  be  re-written 
in  the  desired  form: 


R(r)  =  _  _  x 
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(27) 


where  A(r)  =  2 nr  x  ct/2  is  the  bottom  scattering  area  at 
a  distance  of  r.  9eg(r)  =  (2/n)9eff(r)  =  \jHj  ( nHQ ) 
=  0.637 9eff(r).  Taking  a  logarithm  for  both  sides  of  (27),  we 
have  the  reverberation  level  (in  dB)  as 

RL(r)  =  lOlog l0A(r)+BBS(9eg)  -  2TL(r) 

—  attw  2r  —  6.0,  (28) 


or 


III.  THEORETICAL  MODELS  FOR  THE  BBS{0)  RLir)  =  lOloglo(cT7i2/40//r2)  +  BBS(9eg) 

DERIVATION  FROM  REVERBERATION  DATA  -  attw  x  2 r,  (29) 


In  isovelocity  shallow  water,  the  cycle  distance  of  the 
nth  mode-ray  can  be  expressed  by  Sn  =  2///tan(0„).  Due  to 
mode  stripping,  lower  modes  with  small  grazing  angles 
dominate  the  sound  field  in  the  three-halves  law  region 
where  we  have  tan0„«0„.  Contributions  from  higher- 
order  modes  with  grazing  angles  larger  than  the  critical 
angle  can  be  neglected.  Using  Eq.  (13)  and  Eq.  (10),  the 
sound  intensity  in  the  three-halves  law  region  can  be 
expressed  by: 


Ic{r)tt—e 


Hr 


e~ &Me2d9 
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a  vr 
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for  7-3  >  r  >  7'i . 


Or  \ 

—9c  = - x  29efr(r) 

H  I  Hr  effy  1 


(25) 


Here  9ejf(r)  is  the  effective  angle  of  sound  propagation 
defined  by41 


Mr) 


e-(Qr/H)<?dB 

0 


1  I Hn 
2\j-&- 


for  7-3  >  7-  >  7*2  ,  (26) 


where  the  transition  distances,  74,  r2  and  r2 ,  are  defined  in 
Eq.  (18).  The  angular  spectrum  of  sound  propagation  in  the 
isovelocity  SW  is  a  Gaussian  function.  It  can  be  replaced  by 
a  sharp  cut  at  9eff.  The  sound  waves  propagate  within  an 
effective  grazing  angle  range  of  ±9eff(r);  9ejf(r)  depends  on 
the  water  depth  H,  propagation  distance  r  and  seabed  reflec¬ 
tion  loss  factor  Q. 


where  BBS{9eg )  =  lOlog,^^)"]  =  101og10 [^9eff(r))"] 
is  the  bottom  backscattering  strength  at  9eg(r)  that  equals  to 
0.637 9eff(r).  attw  is  the  sound  absorption  in  water  (dB/m). 
Equation  (28)  is  very  similar  to  the  boundary  reverberation 
formula  suitable  to  short-range  or  deep  ocean  modeling  from 
the  ray  method  (without  considering  waveguide  multipath 
effect).39  Equations  (27)-(29)  show  an  interesting  result:  in 
the  mode  stripping  (three-halves  law)  region,  although  the 
SW  reverberation  is  produced  from  the  bottom  scattering  by 
all  modes  or  rays  with  different  incident  and  scattering  angles, 
the  complex  multi-incident/scattering  angle  interactions  with 
a  bottom  are  equivalently  replaced  by  the  bottom  back- 
scattering  (BBS)  at  a  single  angle  of  9eg  (=0.6377,,//),  i.e.,  the 
reverberation  may  be  expressed  as  a  simple  product  of  the 
transmission  loss  to  and  from  the  insonified  area  and  the 
seabed  backscattering  at  a  single  range-dependent  angle, 
9eg(r).  A  factor  of  1/2  in  Eq.  (27)  for  two-way  propagation 
Ic(r),  or  —6.0  dB  in  Eq.  (28),  accounts  for  the  fact  that  only 
the  down-going  waves,  i.e.,  half  the  energy  of  normal  modes 
in  sound  propagation,  are  considered  as  an  incident  wave 
for  bottom  scattering.  This  factor  is  required  by  the  definitions 
of  the  commonly  used  bottom  scattering  models.  Similarly, 
only  the  up-going  waves  of  the  bottom-scattered  normal- 
modes  are  considered  in  reverberation/scattering  calculations. 

IV.  GEO-ACOUSTIC  MODEL  FOR  SANDY  SEA 
BOTTOMS 

LF  sound  speed  and  attenuation  in  sandy  seabottoms 
have  recently  been  analyzed  and  summarized  from  long- 
range  field  measurements  in  shallow  water.9  Field 
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measurements  conducted  at  20  locations  in  different  coastal 
zones  around  the  world  were  analyzed.  The  sound  attenua¬ 
tions  in  sandy  or  sand-silt  mixture  bottoms,  inverted  from 
different  acoustic  field  characteristics,  exhibit  similar  magni¬ 
tude  and  nonlinear  frequency  dependence  below  2000  Hz  at 
all  of  these  sites  (see  Fig.  2).  The  LF  field-derived  effective 
seabed  geoacoustic  model,  for  both  sound  speed  and  attenua¬ 
tion  in  the  sandy  bottoms,  can  be  well  described  by  the  Biot 
model  as  shown  in  Figs.  2(a)  and  2(b)  by  solid  curves.9  We 
will  use  this  effective  bottom  geoacoustic  model  for  invert¬ 
ing  the  LF&LGA  BBS  from  reverberation  data,  because  half 
of  those  LF  measurements  in  Ref.  9  were  made  at  Chinese 
sea  sites  that  are  the  same  as  or  close  to  the  sites  for  the 
reverberation  measurements  reported  in  this  paper. 

As  Ref.  9  indicates,  in  the  basic  Biot  theory,  three  pa¬ 
rameters  largely  control  the  fluid  motion  that  causes  the  non¬ 
linear  dispersion  that  is  important  in  the  coarser  granular 
sediments.  These  are  the  permeability  ks,  the  porosity  /?,  and 
the  tortuosity  oto  which  defines  the  necessary  amount  of 
“added  mass”  in  the  Biot-Stoll  model.  These  three  parame¬ 
ters  define  the  pore-size  parameter.  Data/model  comparisons 
(trials)  tell  us  that  the  LF  sound  speed  ratio  is  very  sensitive 
to  the  porosity  of  sediment  and  is  less  sensitive  to  the  perme¬ 
ability  and  tortuosity.  The  LF  sound  attenuation  is  very  sen¬ 


sitive  to  the  permeability  of  the  sediment  and  is  less 
sensitive  to  the  porosity  and  tortuosity.  Thus,  the  porosity  of 
0.42  and  the  permeability  of  1.0  x  10-11  were  easily  deter¬ 
mined  as  the  values  that  provide  the  best  match  between  the 
Biot  model  and  the  LF  inverted  attenuations  and  speed  ratios 
shown  in  Fig.  2,  as  well  as  the  average  LF  sound  attenuations 
shown  in  Fig.  1 1  of  Ref.  9.  The  other  eleven  Biot  parameters 
are  same  as  those  from  the  SAX99  measurement.  The  stand¬ 
ard  deviation  of  the  permeability,  inverted  from  the  LF  field 
measurements  is  less  than  0.05  x  10— 1 1 ;  The  porosity  devia¬ 
tion  is  less  than  0.015. 

The  bottom  reflection  factor  Q  in  Eqs.  (13),  (27),  and 
(29)  can  be  expressed  by9,1 1,36-38 


Q(f)  =  0.0366 


(' cj/ch)(ph/pw ) 

,  ,  o]  3/ 

i  -  (cw/chy 


(30) 


where  oq,  is  the  sound  attenuation  (dB/m)  in  the  bottom  and  / 
is  the  frequency  in  kHz.  cw  and  c*  are  the  sound  speed  (km/s) 
in  the  water  and  the  bottom,  respectively.  pw  and  ph  are  the 
density  (g/cm3)  in  the  water  and  the  bottom,  respectively. 

Inputting  the  LF-field  inverted  average  sound  attenua¬ 
tion  and  speed,  shown  in  Figs.  2(a)  and  2(b)  by  solid  lines, 
into  Eq.  (30),  we  have  the  average  bottom  loss  parameter  Q 
vs  frequency  shown  in  Fig.  3  by  solid  line.  The  LF  fit  param¬ 
eters  for  the  Biot  model  are  listed  in  Table  X  of  Ref.  9.  In 
contrast  to  the  Hamilton  model,42’43  which  should  result  in  a 
frequency-independent  Q,  Fig.  3  shows  an  interesting  fre¬ 
quency  dependence  of  the  bottom  reflection  loss  at  LGA:  it 
increases  with  increasing  frequency,  reaches  a  maximum 
around  the  JR  (Jackson  and  Richardson)  characteristic  fre¬ 
quency  of  2.88  kHz  [from  Eq.  (10.52)  of  Ref.  4;  it  is  differ¬ 
ent  from  the  Biot  characteristic  frequency],  then  decreases 
with  frequency.  This  frequency  dependence  of  the  LGA  bot¬ 
tom  reflection  loss  affects  the  two-way  sound  propagation  in 
reverberation  as  well  as  the  frequency  dependence  of  the 
reverberation-derived  bottom  scattering.  Possible  variations 
of  the  Q  values,  caused  by  standard  deviations  of  the  LF 


LF  attenuation  data-Biot  model  comparison,  21  locations 
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FIG.  2.  (Color  online)  Compar  ison  of  the  LF  field-inverted  seabed  acoustic  pa-  FIG.  3.  (Color  online)  The  seabed  reflection  loss  factor  Q  as  a  function  of 

rameters  with  the  Biot  model  (Ref.  9):  (a)  sound  attenuation,  (b)  sound  speed.  frequency  for  sandy  bottoms,  derived  from  the  LF  field  measurements. 
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field-inverted  permeability  and  porosity,  are  plotted  in  Fig.  3 
by  dashed  lines. 

V.  REVERBERATION  MEASUREMENTS  AND 
REVERBERATION-DERIVED  BBS(6) 

A.  Reverberation  measurements 

A  quality  database  of  reverberation  is  absolutely  essen¬ 
tial  in  derivation  of  bottom  scattering  strength  from  reverber¬ 
ation  measurements.  The  “quality”  here  means  that  the 
reverberation  data  are  in  a  broadband,  obtained  from 
carefully  calibrated  measurements.  The  data  have  been 
normalized  to  the  source  level  (SL),  and  have  high 
reverberation-to-noise  ratio.  However,  to  get  wideband  SL- 
normalized  RL  is  a  delicate  task  that  can  be  subject  to  error. 
A  prior  paper  introduced  a  simple  measurement  method 
that  could  avoid  the  signal  overflow  and  saturation  caused 
by  a  powerful  explosive  source.11  The  SL-normalized  RL 
data  in  the  three-halves  law  region  from  this  method  are 
used  in  this  paper  to  derive  the  BBS  (after  the  RL  data  are 
corrected  by  the  Thorp  formula  for  the  sound  absorption  in 
the  ocean  ). 

Reverberation  data  were  collected  using  1000-g  explo¬ 
sive  charges  from  4  sites:  two  from  the  Yellow  Sea  (YS); 
two  from  the  East  China  Sea  (ECS).  During  the  measure¬ 
ments,  the  sea  surfaces  were  relatively  calm.  The  sea  bot¬ 
toms  at  the  three  sites  were  very  flat; 40  the  depth  deviation 
was  less  ±1  m.  The  fourth  site  was  nominally  flat.11'45  The 
average  sound  velocity  profiles  (SVP)  for  the  4  sites,  meas¬ 
ured  at  the  receiving  array  sites,  are  plotted  in  Fig.  4. 

The  SVPs  at  two  sites  (YS  site  I  and  ECS  site  I),  shown 
in  Fig.  4  by  solid  lines,  were  isovelocity.  The  SVP  at  ECS 
site  11  (ASIAEX01),  shown  by  a  dashed  line,  is  quasi¬ 
isovelocity.  At  this  site,  when  explosive  charges  were 
detonated  at  a  depth  of  50  m,  RLs  received  by  different 
hydrophones  located  at  different  depths  did  not  show  any 
depth  dependence,  i.e.,  the  case  was  close  to  a  Pekeris 
model.  The  SVP  at  YS  site  II  (YS’96),  shown  by  a  dotted 
line,  has  a  strong  thermocline.  When  both  sound  source  and 


receiver  are  below  the  thermocline,  the  YS’96  site  can 
approximately  be  treated  as  a  Pekeris  model,  because  numer¬ 
ical  simulations  show  that  the  maximum  RL  difference 
between  a  real  SVP  and  the  Pekeris  model  is  less  than  2  dB 
for  frequencies  larger  than  200  Hz. 


B.  Using  the  Lambert  law  and  the  Biot  geo-acoustic 
model  to  derive  BBS 

The  semi-empirical  Lambert  law  has  widely  been  used 
for  comparisons  with  bottom  scattering  measurements  and 
reverberation  modeling,  although  it  has  no  physical  basis. 
Just  for  the  sake  of  comparison  with  the  commonly  pub¬ 
lished  Lambert  coefficients,46  this  section  will  use  the 
Lambert  law  (n  =  k  +  m  =  2),  reverberation  data  and  the 
Biot  model  to  derive  the  Lambert’s  scattering  coefficient, 
101og10  ///M,  with  respect  to  frequency. 

Figures  5(a)  and  5(b)  show  the  initial  reference  rever¬ 
beration  level  (IRRL)  at  2  and  4  second  after  an  explosive 
source  is  detonated,  respectively. 1 1  The  IRRL  data  were 
obtained  from  4  sites.  Equation  (24)  shows  that  in  the 
three-halves  law  region,  if  the  bottom  scattering  obeys  the 
Lambert  law  («  =  2),  the  reverberation  intensity  is  independ¬ 
ent  of  water  depth,  H.  Using  the  IRRL  data  shown  in  Figs. 
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FIG.  4.  (Color  online)  Sound  velocity  profiles  at  four  sites  for  reverberation 
measurements. 


FIG.  5.  (Color  online)  Reverberation  level  vs  frequency  from  4  sites:  (a)  at 
2  s;  (b)  at  4  s. 
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FIG.  6.  (Color  online)  The  Lambert’s  scattering  coefficient  derived  from 
the  Biot  model  and  RLs  at  (a)  2  s,  (b)  4  s. 


5(a)  and  5(b),  the  Q  values  in  Fig.  3  and  the  closed-form 
expression  of  Eq.  (24),  we  get  the  Lambert  scattering  coeffi¬ 
cients  shown  in  Figs.  6(a)  and  6(b).  The  results  show  that,  if 
the  geoacoustic  model  of  sandy  bottoms  follows  the  Biot 
model,  the  LF  Lambert  scattering  coefficient  exhibits  strong 
frequency  dependence.  Figures  6(a)  and  6(b)  show  that 
below  the  JR  characteristic  frequency,  it  exhibits  frequency- 
cubed  dependence  (ex  /3)  in  the  150-2500  Hz  range.  Aver¬ 
aged  Lambert  scattering  strength  from  Figs.  6(a)  and  6(b) 
can  approximately  be  expressed  by 

lOlogio  Blm  ~  ~  16-4  +  101oglo(//1000)3°(dB) 

for  150-2500  Hz.  (31) 


C.  Using  the  closed  form  expression  for  RL  and  the 
effective  Biot  model  to  derive  BBS(0)  with  arbitrary 
angular  dependence 

If  both  wideband  TL  and  RL  are  well  measured,  or  if 
both  RL  data  and  the  seabed  geoacoustic  model  are  avail¬ 
able,  then  Eqs.  (28)  and  (29)  can  be  used  to  derive  an  angular 


dependence  of  the  LF  BBS (0).  In  this  section,  the  wideband 
RL  data,  obtained  from  three  sites  with  sandy  bottoms,  and 
the  average  Q  values  shown  in  Fig.  3,  are  used  in  Eq.  (29)  to 
derive  the  BBS  vs  angle  and  frequency.  The  RL  data  from 
the  YS  site  II  are  not  used,  because  the  RL  data  at  this  site 
exhibit  an  uncommon/unknown  water  volume  scattering  at 
long-range  (around  20  s). 

Equation  (18)  shows  that  the  three-halves  law  region  is 
within  two  transition  ranges,  r 2  and  rj,  .  That  is,  the  RL  data 
used  to  derive  BBS(d)  in  this  section  should  satisfy  two  con¬ 
ditions:  (1)  A  high  reverberation/noise  ratio  (>6  dB);  (2) 
Reverberation  time  t  between  2/2 /c  and  2/q/c  (c  is  the  sound 
speed  in  the  water  column).  Based  on  Eqs.  (26)  and  (29), 
r 2,  /'3,  H  (water  depth)  and  Q  (the  bottom  loss  factor)  define 
a  possible  angular  range  of  the  reverberation-derived  BBS (9) 
for  a  given  experimental  site. 

( 1 )  For  the  YS  site  I,  where  the  water  depth  is  28.5  ±  1 .0  m,  the 

RL-derived  BBS (6)  as  a  function  of  grazing  angle  for 
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FIG.  7.  (Color  online)  Angular  and  frequency  dependence  of  the  reverbera¬ 
tion-derived  BBS  from  the  YS  Site  I.  (a).  BBS  vs  angle  for  eight  frequen¬ 
cies,  (b)  BBS  vs  frequency  at  4°,  6°,  8°,  and  10°. 
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different  frequencies  is  plotted  in  Fig.  7(a)  by  dashed  lines. 
Figure  7(a)  covers  eight  center  frequencies  with  1/3  octave 
bands  (from  the  bottom  line  in  the  figure):  315,  500,  800, 
1000,  1250,  1600,  2000,  and  2500  Hz"  For  highlighting  the 
BBS  characteristics  at  LGA,  the  horizontal  axis  of  Fig.  7(a) 
is  in  a  logarithmic  scale.  Data  fits  show  that  the  BBS  at  very 
low  grazing  angle  (l°—4°)  in  the  800-2500  H  range  is  pro¬ 
portional  to  ks  sin4  0,  and  ks  «  7  —  70.  A  straight  dashed  line 
of  13.6  +101og(sin4  9)  dB  is  shown  in  Fig.  7(a);  it  corre¬ 
sponds  to  ks  =  23.  For  comparison,  the  semi-empirical 
Lambert  law,  —27  +  10  logfsirr  6),  is  plotted  by  a  straight 
solid  line. 

The  values  of  BBS(9)  at  4°,  6°,  8°  and  10°  for  differ¬ 
ent  frequencies  at  the  YS  Site  I  are  shown  in  Fig.  7(b). 
Between  315  Hz  and  2500  Hz,  this  can  be  approximately 
expressed  by 

BBS{6)  «  -38.8  +  101oglo(/yi000)2-2  (dB).  (32) 


(2)  For  the  ECS  site  I,  where  the  water  depth  is  40.0 ±  1.0m, 
the  RL-derived  BBS (9)  as  a  grazing  angle  for  different 
frequencies  is  plotted  in  Fig.  8(a)  by  dashed  lines.  Figure 
8(a)  covers  ten  center  frequencies  with  1/3  octave  bands 
(from  the  bottom  line  in  the  figure):  315,  400,  500,  630, 
800,  1000,  1250,  1600,  2000,  and  2500  Hz.  The  BBS  at 
2000  Hz  for  very  low  grazing  angle  (1°— 4°)  is  about 
100  sin4  9  shown  by  a  straight  dashed  line.  All  of  the 
BBS(9)  at  4°,  6°,  8°  and  10°  for  different  frequencies 
from  the  ECS  site  I  are  shown  in  Fig.  8(b).  Between  315 
Hz  and  2500  Hz,  it  can  be  approximately  expressed  by 

BBS{9)  «  -34.5  +  101og10(/'/1000)2'°  (dB).  (33) 

(3)  For  the  ECS  site  11  (ASIAEX01),  where  the  seabed  is 
nominally  flat,11’45  with  an  average  water  depth  of  104  m 
over  the  extent  of  the  reverberation  measurements,  the 
RL-derived  BBS (9)  as  a  grazing  angle  for  different 
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FIG.  8.  (Color  online)  Angular  and  frequency  dependence  of  the 
reverberation-derived  BBS  from  the  ECS  Site  I.  (a).  BBS  vs  angle  for  ten 
frequencies,  (b)  BBS  vs  frequency  at  4°,  6°,  8°,  and  10°. 


FIG.  9.  (Color  online)  Angular  and  frequency  dependence  of  the 
reverberation-derived  BBS  from  the  ASIAEX01  Site  I.  (a).  BBS  vs  angle 
for  fourteen  frequencies,  (b)  BBS  vs  frequency  at  4°,  6°,  8°,  and  10°. 
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frequencies  is  plotted  in  Fig.  9(a)  by  dashed  lines.  Figure 
9(a)  covers  fourteen  center  frequencies  with  1/3  octave 
bands  (from  the  bottom  line  in  the  figure):  200,  300,  400, 
500,  600,  700,  800,  1000,  1250,  1500,  1750,  2000,  2250, 
and  2500  Hz.  The  scattering  coefficient  at  2000  Hz  for 
low  grazing  angle  ( 1° — 4°)  is  about  4  sin3  9  as  shown  by  a 
straight  dashed  line.  All  BBS(0)  at  4°,  6°,  8°  and  10°  for 
different  frequency  from  the  ECS  site  II  are  shown  in 
Fig.  9(b).  Between  200  Hz  and  2500  Hz,  it  can  approxi¬ 
mately  expressed  by 

BBS{6)  «  -34.7  +  101og10(/'/1000)2'7(dB).  (34) 


VI.  SUMMATION  AND  DISCUSSION 

The  angular  spectrum  (energy  flux)  method  for  calculat¬ 
ing  SW  reverberation  has  briefly  been  introduced  in  this  pa¬ 
per.  It  is  based  on  the  WKB  approximation  to  the  normal¬ 
mode  solution  for  wave  equations  (ray-mode  analogies).  In 
this  method,  two-way  sound  propagation  and  bottom  scatter¬ 
ing  are  treated  as  weighting  functions  in  the  angular  domain. 
As  an  alternative  method  to  SW  reverberation  modeling,  it  is 
simple  and  intuitive. 

The  closed-form  expressions  for  SW  reverberation  in  the 
isovelocity  waveguide,  derived  from  the  energy  flux  method, 
show  that  the  complex  seabed  scattering  by  all  modes/rays  with 
different  incident  and  scattering  angles  can  equivalently  be  rep¬ 
resented  by  seabed  back-scattering  at  a  single  angle  (that 
depends  on  range,  seabed  reflection  loss  and  water  depth).  This 
equivalent  relationship  offers  a  simple  method  to  derive 
LF&LGA  BBS  from  reverberation  measurements  [see  Eqs. 
(26) — (29)] .  The  method  is  valuable  because  it  allows  seabed 
backscattering  strength  to  be  inverted  at  frequencies  and  angles 
where  no  practical  method  of  measurements  exists. 

The  equivalent  relationship  and  the  LF  field-derived 
seabed  geo-acoustic  model  (Biot  parameters)  are  used  to 
extract  the  bottom  backscattering  strength  in  a  frequency 
band  of  200-2500  Hz  and  in  a  grazing  angle  range  of 
1.1°-14.0°  from  SW  reverberation  measurements  at  three 
sites  with  sandy  bottoms.  The  results  are  shown  in  Figs.  7(a), 
8(a),  and  9(a).  From  these  figures,  many  observations  can  be 
made:  (1)  In  general,  in  the  500-2500  Hz  range,  the  SW 
reverberation-derived  BBSs  at  low  grazing  angles  are  larger 
than  the  values  predicted  by  the  semi-empirical  Lambert  law. 
(2)  For  grazing  angles  that  are  larger  than  about  4°,  the  BBS 
has  weaker  angle  dependence  than  the  Lambert  law.  (3)  For 
the  high  end  of  frequencies  analyzed  in  this  paper,  if  grazing 
angles  are  smaller  than  around  4°,  the  BBS  decreases  faster 
than  the  Lambert  law  with  decreasing  angles,  the  angular 
exponents  may  reach  3^1.  (4)  In  a  frequency  range  of 
200-2500  Hz  (below  the  JR  characteristic  frequency),  the 
LF&LGA  BBS  exhibits  strong  frequency  dependence  for  all 
three  sites.  At  low  grazing  angles  of  4°-10°  it  can  approxi¬ 
mately  be  expressed  by  BBS  =  —  A  +  101og10(//1000)B 
(dB),  where  A  is  a  range  of  34.5-38.8,  B  is  in  a  range  of 
2. 0-2. 7. 

The  characteristics  of  the  reverberation-derived  LF&LGA 
bottom  backscattering  strength  need  a  physical  explanation. 
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The  results  might  be  used  to  compare  with  the  well-developed 
HF  seabed  scattering  models  summarized  by  Jackson  and 
Richardson,6  and  to  see  if  those  HF  seabed  scattering  mod¬ 
els  are  applicable  at  low  frequencies.  It  also  might  be  used 
to  analyze  the  physical  mechanisms  of  bottom  scattering, 
including  possible  estimates  of  the  seabed  roughness  spec¬ 
tra  or  sediment  inhomogeneities  at  low  frequencies.47 

The  theoretical  expressions  used  to  invert  the  BBS,  Eqs. 
(24),  (26)-(29),  are  limited  to  the  three-halves  law  region 
where  r?,  >  r  >  r2-  In  addition  to  this  limitation,  it  is  often 
hard  to  get  quality  reverberation  data  at  long  range  because 
of  ocean  background  noise  interference;  particularly  at  lower 
frequencies.  These  two  factors  limit  the  possible  angular 
ranges  of  the  SW  reverberation-derived  BBS  in  this  paper. 

The  LF  field-inverted  seabed  physical  parameters  for 
the  Biot  model910  and  the  SW  reverberation-derived 
LF&LGA  seabed  scattering  in  this  paper  are  limited  to  sandy 
or  sand-silt  mixture  seabottoms.  Both  seabed  reflection  and 
scattering  models  for  the  BBS  derivations  are  described  by 
the  approximate  expressions,  Eqs.  (13)  and  (19).  The  shear 
wave  effects  are  assumed  negligible  in  the  frequency  range 
for  the  BBS  derivations.  Thus,  the  reverberation-derived 
LF&LGA  BBSs  in  this  paper  are  approximation  values. 

Equations  (26)-(29),  used  to  derive  the  BBS  in  this  pa¬ 
per,  show  how  a  variation  of  the  bottom  loss  factor  Q  would 
quantitatively  affect  the  values  of  the  effective  angle  and 
reverberation-derived  BBS.  The  Q  values  in  Fig.  3  obtained 
from  the  Biot  model  exhibit  very  interesting  frequency  de¬ 
pendence  at  low  grazing  angles.  The  reflection  loss  amplitude 
and  peak  value  frequency  (around  the  JR  characteristic  fre¬ 
quency)  are  a  complex  function  of  the  Biot  parameters,  and 
are  particularly  sensitive  to  the  permeability  ks,  the  porosity 
(S,  and  the  tortuosity  ao-  The  uncertainty  of  the  reverberation- 
derived  BBS  in  this  paper  will  depend  on  the  reliability  of  the 
effective  Biot  model  (average  parameters)  inverted  from  LF 
field  measurements.9,10  Our  preliminary  data-model  compari¬ 
sons  also  show  that  around  (or  higher  than)  the  JR  character¬ 
istic  frequency,  the  reverberation-derived  LGA  BBS  exhibits 
much  weaker  frequency  dependence  than  LF  BBS,  and  it 
sometimes  decreases  with  increasing  frequency.  Thus,  more 
work  needs  to  be  done  in  the  future  to  investigate  how  the 
Biot  parameters  and  their  uncertainties  affect  the  seabed 
backscattering  strength  and  its  frequency  dependence,  includ¬ 
ing  possible  effects  of  shear  wave. 

In  order  to  compare  with  the  inversion  method  used  in 
this  paper,  one  may  simultaneously  measure  both  TL  and  RL 
from  same  sites  and  use  Eq.  (28)  to  derive  the  LF&LGA 
BBS.  It  will  be  desirable  in  the  future  to  directly  integrate 
physics-based  seabed  reflection/scattering  models4  into  the 
energy  flux  method  for  SW  reverberation  modeling,  and  to 
do  data/model  comparisons,  then  from  reverberation  data  to 
invert  the  physical  parameters  of  seabed  scattering  as  well  as 
LF&LGA  BBS  vs  frequency  and  angle.43 
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Several  physics-based  seabed  geoacoustic  models  (including  the  Biot  theory)  predict  that  compres- 
sional  wave  attenuation  in  sandy  marine  sediments  approximately  follows  quadratic  frequency  de¬ 
pendence  at  low  frequencies,  i.e.,  12  ~  kfn  (dB/m),  n  =  2.  A  recent  paper  on  broadband  geoacoustic 
inversions  from  low  frequency  (LF)  field  measurements,  made  at  20  locations  around  the  world,  has 
indicated  that  the  frequency  exponent  of  the  effective  sound  attenuation  n  ss  1.80  in  a  frequency  band 
of  50-1000  Hz  [Zhou  et  al.,  J.  Acoust.  Soc.  Am.  125,  2847-2866  (2009)].  Carey  and  Pierce  hypothe¬ 
size  that  the  discrepancy  is  due  to  the  inversion  models’  neglect  of  shear  wave  effects  [J.  Acoust.  Soc. 

Am.  124,  EL271-EL277  (2008)].  The  broadband  geoacoustic  inversions  assume  that  the  seabottom  is 
an  equivalent  fluid  and  sound  waves  interact  with  the  bottom  at  small  grazing  angles.  The  shear  wave 
velocity  and  attenuation  in  the  upper  layer  of  ocean  bottoms  are  estimated  from  the  LF  field-inverted 
effective  bottom  attenuations  using  a  near-grazing  bottom  reflection  expression  for  the  equivalent 
fluid  model,  derived  by  Zhang  and  Tindle  [J.  Acoust.  Soc.  Am.  98,  3391-3396  (1995)].  The  resultant 
shear  wave  velocity  and  attenuation  are  consistent  with  the  SAX99  measurement  at  25  Hz  and 
1000  Hz.  The  results  are  helpful  for  the  analysis  of  shear  wave  effects  on  long-range  sound  propaga¬ 
tion  in  shallow  water.  ©2012  Acoustical  Society  of  America. 

[http://dx.doi.org/10.1121/L4765078] 
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I.  INTRODUCTION 

Low  frequency  (LF)  (compressional)  sound  velocity  and 
effective  attenuation  in  sandy  seabottoms  have  recently  been 
analyzed  and  summarized  from  long-range  field  measure¬ 
ments  in  shallow  water  (SW).1  The  measurements  were  con¬ 
ducted  by  different  investigators  at  20  locations  in  different 
coastal  zones  around  the  world.  The  effective  seabed  attenua¬ 
tions,  inverted  from  these  acoustic  field  measurements  with 
different  characteristics,  exhibit  similar  magnitude  and  nonlin¬ 
ear  frequency  dependence  below  2000  Hz  at  all  of  these  sites, 
as  shown  in  Fig.  1.  The  average  LF-fie ld-inverted  sound 
attenuation  values  from  sites  1-20  are  listed  in  Table  I  as  a 
function  of  frequency.  The  numbers  of  the  available  data  sets 
for  the  average  are  also  listed  in  Table  I.  Using  power  law  fit¬ 
ting,  the  average  LF  effective  attenuation  &2E  in  sandy  and 
sand-silt  mixture  seabottoms  can  approximately  be  expressed 
as  the  following  nonlinear  frequency  dependence: 

a2£  =  (0.3 67 ±0.011) (//l 000) ( 1 80 1  * 0020)  dB/m, 

for  50- 1000  Hz,  (1) 

»2£  = (0 . 3 72  ±  0 .022)  (// 1 000)  (I-809±0  032)  dB/m, 

for  50 -500  Hz,  (2) 

where /is  frequency  in  units  of  Hz.  The  corresponding  aver¬ 
age  sound  velocity  ratio  at  the  bottom-water  interface  in  the 
50-600  Hz  range  is  1.061  ±  0.009.1 


a)Also  at:  State  Key  Laboratory  of  Acoustics,  Institute  of  Acoustics,  Chinese 
Academy  of  Sciences,  Beijing  100190,  China. 


Numerical  simulations2  show  that  that  both  LF  field- 
derived  sound  velocity  and  attenuation  in  sandy  seabottoms 
can  closely  be  predicted  by  the  Biot-Stoll  model,3-14  the 
Buckingham  VGS  (viscous-grain-shearing)  model,15-17  the 
Chotiros-Isakson  BICSQS  (Biot-Stoll  with  squirt  flow  and 
shear)  model,18  and  the  Pierce-Carey/Williams  simpler 
model19'20  when  the  physical  parameters  of  these  seabed 
acoustic  models  are  known.  However,  there  remain  impor¬ 
tant  unanswered  questions.  One  of  these  involves  the  possi¬ 
ble  shear  wave  effects  that  were  assumed  to  be  unimportant 
in  those  LF  measurements  and  inversion  methods.  Another 
one  is,  why  is  the  frequency  exponent  of  sound  attenuation 
in  Eqs.  (1)  and  (2)  around  1.8,  not  close  to  2  that  should  be 
“a  direct  consequence  of  the  Biot  theory”?17-21  There  might 
be  different  answers  to  the  second  question.  Carey  et  al. 
hypothesized  that  “The  discrepancy  can  be  explained 
because  the  inverse  analysis  inferences  were  made  with  the 
neglect  of  an  additional  attenuation  mechanism  where  gen¬ 
erated  lower  shear  waves  carry  energy  downwards  out  of 
the  waveguide.”21  Following  this  hypothesis,  this  paper 
attempts  to  use  the  LF  field-inverted  effective  sound  attenu¬ 
ation  shown  in  Table  I  to  estimate  a  possible  range  of  the 
average  shear  wave  velocity  and  attenuation  in  the  upper 
layer  of  sandy  bottoms. 

II.  ESTIMATION  OF  SHEAR  WAVE  VELOCITY  AND 
ATTENUATION  FROM  THE  LF  FIELD  MEASUREMENTS 

The  LF  field  geo-acoustic  inversions  in  Ref.  1  assumed 
that  sea  bottoms  were  a  half-space  equivalent  fluid  model; 
i.e.,  the  shear  wave  was  negligible.  Thus,  inverted  sound 
attenuations  (and  related  Biot  parameters)  in  the  bottoms 
should  be  considered  as  equivalent/effective  values,  ot2E(f)- 
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FIG.  1.  (Color  online)  Effective  LF  sound  attenuations  in  sandy  seabottoms 
(Ref.  1).  All  data  symbols  and  a  dotted  line  are  from  the  LF  field  measure¬ 
ments;  a  solid  line  is  obtained  from  the  Biot-Stoll  geoacoustic  model. 


A.  Shear  wave  effects  on  the  near-grazing  acoustic 
reflection  and  the  equivalent  sound  attenuation 


The  LF  sound  velocity  and  attenuation  in  Ref.  1  were 
inverted  from  the  long-range  sound  field  where  lower-order 
normal  modes  interact  with  an  upper  layer  of  the  bottom  at 
low  grazing  angles  (around  a  wave  length  thick).  According 
to  the  Wentzel-Kramers-Brillouin  approximation,  the  nth 
modal  attenuation  coefficient  in  a  SW  waveguide,  //,,  can  be 
expressed  approximately  in  terms  of  the  bottom  reflection 
coefficient  V(9), 


nmol 

Sn{9n)  +  <5  n 


(3) 


where  Sn  is  the  cycle  distance  of  the  nth  mode-ray,  <5„  is  a 
beam  displacement  on  the  bottom  reflection,  9n  is  the  graz¬ 
ing  angle  of  the  nth  mode  at  the  water-bottom  interface.  If 
the  water  depth  H  and  wave  number  k  satisfy  kH  1 ,  <5„  is 
negligible.  The  plane  wave  reflection  coefficient  from  the 
seabed  (fluid  or  elastic)  can  always  be  expressed  by  its  mod¬ 
ulus  and  phase 

V(9)  =  |V(0)|exp[/<K0)]-  (4) 

At  low  grazing  angle, 22-26 


imi  * 

«  1  -Q9^e-Qe  (0<  1), 

(5) 

-n  +  P9  (0<  1). 

(6) 

The  quantity  Q  governs  the  modal  attenuation;  it  is 
related  to  the  bottom  reflection  loss  defined  by  2Olog|V'(0)| 
=  8.686 Q9  in  dB/radian;  P  is  a  near-grazing  reflection  phase 
shift  that  is  introduced  by  an  imaginary  ideal  pressure- 
release  boundary  located  a  distance  of  P/4n  wavelength 
below  the  true  seabed  boundary,  i.e.,  a  “hidden  depth” 
AH  =  (P/4n)L  The  modal  decay  factors  jln  control  the 
sound  propagation  loss  in  SW  waveguides.  Based  on 
Eqs.  (3)  and  (5),  the  shear  wave  must  affect  long-range  SW 
propagation  through  the  near-grazing  bottom  reflection  fac¬ 
tor  Q  and  P. 

The  characteristics  of  the  seabed  reflection  loss  at  low 
grazing  angles  were  used  by  Tindle,  Zhang,  and  Chapman  to 
treat  the  shear  wave  effects25-29  that  will  be  briefly  described 
next. 

The  seabed  is  assumed  as  a  homogeneous  poro-elastic 
half-space  with  density  p2,  compressional  wave  velocity  C2, 
and  shear  wave  velocity  cs.  The  water  layer  above  the  seabed 
has  a  density  p1  and  sound  velocity  c \  (c2  >  c\,  cs  <  C\). 
Attenuation  of  the  compressional  waves  and  shear  waves 
in  the  seabed  is  introduced  by  allowing  their  velocities 
to  become  complex,  c2  — >  c2(  1  —  ie 2)  and  cs  — >  Cj(l  —  ies ) 
with  £2  and  es  small  dimensionless  numbers,  i.e.,  0.2 
=  (2nf/c2)£2  and  a s  =  (2nf/cs)es  in  Np/m. 

The  main  idea  introduced  by  Tindle,  Zhang,  and  Chap¬ 
man  is  to  have  a  fluid  seabed  model  that  is  equivalent  to  an 
elastic  seabed  model,  to  let  three  equivalent  parameters 
(p2£)  C2e,  and  o.2e)  do  the  work  of  five  {p2,  C2,  «2,  cs,  and 
0£j),  and  to  have  the  same  bottom  reflection  parameters  (Q  and 
P)  and  the  same  nonnal-mode  decay  factors  /?„  in  a  shallow- 
water  waveguide.  For  a  low  shear  velocity  ocean  bottom  with 
cs  <  c„.  es  <  1,  £p  <  1,  ignoring  £2,  £2,  epc ],  esc5  and 
higher  order,  and  leaving  the  critical  angle  (equivalent  C2) 
unchanged 


C2E  =  c2.  (7) 

Tindle,  Zhang,  and  Chapman  derived  the  following  approxi¬ 
mation  expressions  for  the  equivalent  sound  attenuation  («2e), 
density  p2E,  near-grazing  bottom  reflection  factors  (Q),  phase 
shift  (P),  and  the  hidden  depth  in  terms  of  real  acoustic  param¬ 
eters  and  the  angular  frequency  cu  {=2nf)25~19 


TABLE  I.  The  average  effective  sound  attenuation  from  the  LF  field  measurements. 


/(Hz) 

50 

80 

100 

125 

160 

200 

250 

300 

315 

400 

500 

«2  E 

(dB/m) 

0.00175 

0.00330 

0.00598 

0.00878 

0.01510 

0.01978 

0.02950 

0.04271 

0.04368 

0.07012 

0.10926 

Data  sets 

(7) 

(7) 

(15) 

(5) 

(10) 

(19) 

(9) 

(15) 

(8) 

(19) 

(19) 

/(Hz) 

600 

630 

700 

750 

800 

900 

1000 

1200 

1500 

1600 

2000 

X2E 

(dB/m) 

0.14515 

0.18479 

0.18888 

0.23650 

0.22365 

0.27511 

0.36504 

0.47193 

0.75464 

0.67210 

1.17768 

Data  sets 

(12) 

(9) 

(12) 

(6) 

(12) 

(7) 

(13) 

(8) 

(5) 

(3) 

(4) 
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a2£  —  «2  + 


Acoc2c]{l-c\lclfl\l-c]/c\)ll2 
c5(l-2c2/c2)2 

8a.sc2c3(l  -  c\lc\) 

c\{\-2c2/c\)  ’ 


(8) 


P2E  =  PlU  -Z^s/cl}2 


t  {*2eYc\{\  +  2c\/cI) 
2m2{\  -c\lc22)2 


1  “I 


(9) 


2  =  2 


P2  fj_C2 

Pi  [271/ 

<ty 


a2[l-2(cJ/ci)“] 


2i2 


(g/c2)2 


[i-(Cl/C2)2]3/2 


1  8c, 


+4(^)  \/l~(cs/ci)2  +  --^«s(cs/c i)2 


[1  — 2(cj/ci)2 

[1-(C!/C2)2] 


2l  1/2  ’ 


p  =  2y0  =  2(p2/Pl) 


p 

A H=  L 

471 


[l-2(c,/d)2]2 

[1-(ci/c2)2]1/2' 


(10) 

(11) 

(12) 


In  these  equations  a2  and  oq  are  in  Np/m.  Based  on  these 
equations.  Chapman  clearly  indicates  “In  propagation  condi¬ 
tion,  where  only  Q  and  P  determine  the  held,  there  are  many 
different  seabottoms  (many  families  of  five  parameters)  that 
deliver  the  same  acoustic  effect.  Conversely,  long-range 
inversion  experiments  that  effectively  only  measure  Q  and  P 
cannot  hope  to  produce  unambiguous  values  of  five  seabed 
parameters.  However,  the  above  equations  strongly  suggest 
how  estimated  parameter  values  might  be  correlated.  For 
example,  an  inversion  algorithm  that  underestimates  shear 
velocity — or  perhaps  assumes  zero  shear  velocity — would 
naturally  also  underestimate  the  density  and  overestimate  the 
compressional-wave-attenuation,  to  compensate.  Of  course, 
the  ambiguity  of  these  inversions  could  be  resolved  by  sam¬ 
pling  the  held  at  short  range,  where  the  acoustic  held  con¬ 
tains  information  about  seabed  interactions  at  steeper 
angles.”25  Unfortunately,  the  acoustic  held  at  short  range  is 
often  too  insensitive  to  infer  compression  wave  attenuation 
in  the  seabed;  on  other  hand,  the  inverted  seabed  acoustic  pa¬ 
rameters  that  are  obtained  at  short  range  might  not  represent 
the  average  geoacoustic  characteristics  of  the  entire  large 
area  used  in  the  sonar  performance  modeling.  Thus,  using 
inversions  to  estimate  the  parameters  of  an  effective  seabed 
model  “may  turn  out  to  be  more  useful  for  sonar  perform¬ 
ance  prediction  models.”25 

Based  on  Eq.  (8),  all  the  LF  held  inversions  in  Ref.  1 
that  ignored  the  shear  wave  effect  result  in  an  equivalent  bot¬ 
tom  attenuation,  a2£,  that  is  larger  than  the  real  compres- 
sional-wave-attenuation,  a2.  Tindle,  Zhang,  and  Chapman 
numerically  show  that  the  equivalent  fluid  approximations 
are  capable  of  providing  a  good  approximation  to  the  reflec¬ 
tion  coefficient  of  a  solid  seabed  (with  a  low  shear  wave  ve¬ 
locity),  as  well  as  a  good  approximation  to  SW  long-range 

25—29 

propagation.^ 


B.  Shear  velocity  and  attenuation  in  sandy  and 
sand-silt  mixture  seabottoms,  estimated  from  the  LF 
field-inverted  equivalent  sound  attenuation 

The  derivation  of  shear  velocity  and  attenuation  in  this 
section  is  based  on  the  following  four  conditions  (models, 
measurements,  or  assumptions):  (i)  the  LF  compressional 
wave  attenuation  in  ocean  bottoms  should  follow  quadratic 
frequency  dependence  (ex/20)17  21;  (ii)  the  averaged  LF 
held-inverted  equivalent  sound  attenuation,  a 2e,  and  veloc¬ 
ity,  c2e,  as  a  function  of  frequency  have  been  obtained;  (iii) 
three  equivalent  parameters  (p2£,  c2£,  and  a2^)  in  an  equiva¬ 
lent  fluid  seabed  model  do  the  work  of  five  (p2,  c2,  a2,  cs, 
and  a*)  in  a  poro-elastic  seabed  model  for  the  near-grazing 
bottom  reflection  (consequently  for  SW  long-range  sound 
propagation)  described  by  Eqs.  (7) — (1 1);25  29  and  (iv)  shear 
wave  velocity  and  attenuation  in  sandy/silt  bottoms  obey  the 
Biot-Stoll  geoacoustic  model. 

Numerical  simulations  show  that  when  both  frame  bulk 
and  shear  moduli  are  negligibly  small,  or  two  moduli  are  real 
(have  no  an  imaginary  term),  the  Biot-Stoll  model  does  pre¬ 
dict  a  quadratic  frequency  dependence  at  low  frequencies. 

As  we  know,1  the  higher  the  frequency,  the  larger  the 
uncertainty  of  the  geoacoustic  inversion  in  shallow  water 
would  be.  Yet  at  low  frequencies,  the  sound  attenuation, 
derived  from  several  physics-based  geoacoustic  models, 
exhibits  quadratic  frequency  dependence.  Taking  these  two 
considerations  into  account,  we  will  only  use  the  LF  field- 
inverted  attenuation  data,  a2£</),  between  50  Hz  and  500  Hz 
(that  are  listed  in  Table  I). 

In  the  fall  of  1999,  a  comprehensive  field  experiment 
named  “SAX99”  (sediment  acoustics  experiments)  was  con¬ 
ducted  in  the  Gulf  of  Mexico  (referred  to  as  SAX99 
site).1314  During  the  SAX99  experiment,  the  Biot  parame¬ 
ters  of  the  sediments  were  extensively  measured  by  using 
both  traditional  and  newly  developed  methods.  The  sound 
velocity  was  measured  over  a  frequency  range  of  125  Hz- 
400  kHz,  and  attenuation  was  measured  over  a  frequency 
range  of  2.6-400  kHz  by  scientists  from  different  institu¬ 
tions.  The  data  are  unique  both  for  the  frequency  range 
spanned  at  a  common  location  and  for  the  extensive  environ¬ 
mental  characterization.  All  the  results  on  compressional 
and  shear  velocities  and  attenuations  were  summarized  by 
Williams  el  al. 1 3  Reference  1  shows  that  the  LF  sound  veloc¬ 
ity  and  attenuation  in  sandy  bottoms,  inverted  from  20  sites, 
smoothly  join  with  the  SAX99  benchmark  data  set  in  the 
mid  to  high  frequency  ranges.  Thus,  the  shear  velocities  at 
two  frequencies  of  1000  Hz  and  25  Hz  from  the  SAX9913'14 
will  be  used  in  this  paper  as  reference  values.  Shear  wave 
velocities  were  obtained  from  the  SAX99  experiments  using 
time-of-flight  measurements.1.  At  1000  Hz,  the  average 
shear  velocity  was  122  m/s  with  a  range  of  97  to  147  m/s;  the 
attenuation  of  shear  waves  was  found  to  have  a  mean  of 
30.5  dB/m  with  a  range  of  21  to  40  dB/m  at  1  kHz.  At  25  Hz, 
the  mean  shear  wave  velocity  was  150  m/s  (in  the  upper 
layer  of  about  lm  thick,  it  was  about  lOOm/s;  beneath  this 
layer,  the  maximum  velocity  may  approach  about  200  m/s); 
the  mean  shear  wave  attenuation  was  0.33  dB/m  (in  the 
upper  layer  it  was  0.47  dB/m,  in  the  lower  stratum  it  is 
0.19  dB/m).14 
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In  the  basic  Biot-Stoll  theory,  three  parameters  largely 
control  the  fluid  motion  that  causes  the  nonlinear  dispersion 
that  is  important  in  the  coarser  granular  sediments.  These  are 
the  permeability  ks,  the  porosity  f>,  and  the  tortuosity  ao» 
which  define  the  necessary  amount  of  “added  mass”  in  the 
Biot-Stoll  model.  These  three  parameters  define  the  pore- 
size  parameter.  Data/model  comparisons  (trials)  tell  us  that 
the  LF  sound  velocity  ratio  is  very  sensitive  to  the  porosity 
of  sediment  and  is  less  sensitive  to  its  permeability  and  tortu¬ 
osity.  The  LF  sound  attenuation  is  very  sensitive  to  the  per¬ 
meability  of  the  sediment  and  is  less  sensitive  to  its  porosity 
and  tortuosity.  Thus,  the  porosity  of  0.420  ±  0.015  and  the 
permeability  of  (1.00  ±0.05)  x  10  11  were  easily  deter¬ 
mined  as  the  values  that  provide  the  best  match  between  the 
Biot-Stoll  model  and  the  LF  inverted  attenuations  and  veloc¬ 
ity  ratios.  The  other  eleven  Biot  parameters  were  taken  same 
as  those  from  the  SAX99  measurement.’1  All  the  effective 
Biot  parameters  that  were  used  to  match  the  LF  field- 
inverted  sound  velocity  and  attenuation  are  listed  in  Ref.  1 
(Table  X,  in  the  column  labeled  “LF  fit”).  The  frame  shear 
wave  modulus  /t0  and  shear  log  decrements  Ss  are  not  impor¬ 
tant  for  the  compressional  wave  data-model  comparison  in 
Ref.  1;  however,  they  play  significant  roles  in  determining 
the  shear  wave  velocity  and  attenuation  in  this  paper. 

The  procedures  for  deriving  shear  wave  velocity  and 
attenuation  from  LF  field  measurements  are  as  follows:  (i) 
by  adjusting  a  pair  of  frame  shear  moduli  /(0  and  3S  and 
using  the  Biot-Stoll  model,  calculate  shear  wave  velocity 
and  attenuation  vs  frequency,  i.e.,  cs(f)  and  xs(f)\  (ii)  input 
the  resultant  cs(f)  and  is(f)  into  Eq.  (8)  and  use  the  LF-field 
inverted  equivalent  sound  attenuation  and  velocity,  «2 e(J) 
and  C2{=C2e),  respectively,  to  calculate  a  set  of  real  com¬ 
pressional  sound  attenuation  as  a  function  of  frequency,  i.e., 
0£2 (/);  (iii)  when  012(f)  is  proportional  to/2000,  a  pair  of  /t0 
and  ds,  and  a  pair  of  cs(f)  and  1  s(f)  are  inverted. 

The  third  term  of  Eq.  (8)  indicates  that  there  is  shear 
velocity-attenuation  coupling;  a  decrease  in  shear  velocity 
can  be  compensated  by  an  increase  in  shear  wave  attenua¬ 
tion;  that  is,  there  may  be  many  combinations  of  cs  and  as 
that  deliver  the  same  effect.  Thus,  we  use  the  shear  velocity 
at  25  Hz  or  at  1000  Hz  obtained  from  the  comprehensive 
SAX99  measurements  as  a  constraining  condition.  Based  on 
“historical  bounds”  of  frame  shear  modulus  /.<0  and  shear  log 
decrements  5S  for  sandy/silt  bottoms,30  and  //0  was  varied 
between  2.61  x  107  —  11.9  x  107  at  0.01  x  107  intervals 
with  Ss  between  0-0.15  at  0.0001  intervals,  we  calculate 
shear  velocity  and  attenuation  as  inputs  for  Eq.  (8);  the  LF 
field-inverted  equivalent  bottom  attenuation  0.2 e  is  numeri¬ 
cally  modified  toward  the  compressional  wave  attenuation 
0C2-  The  pair  of  Ho  and  5S  values,  and  the  corresponding  pair 
of  cs(f)  and  as(f)  values,  for  which  the  resultant  sound 
attenuation  «2  has  the  quadratic  frequency  dependence 
(oc/2  000),  are  both  readily  uniquely  determined.  The  result¬ 
ant  cs  and  a,s  as  a  function  of  frequency  are  shown  in  Figs. 
2(a)  and  2(b).  Solid  lines  (case  A)  in  these  two  figures  are 
obtained  by  using  cs  =  122  m/s  as  a  constraining  condition. 
It  was  obtained  by  the  NRL  team  from  the  SAX99  site  at 
1 000  Hz. 1  The  dashed  lines  (case  B)  are  obtained  by  using 
cs=  150  m/s  as  a  constraining  condition.  It  was  obtained  at 
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FIG.  2.  (Color  online)  LF  field  inverted  shear  wave  property  in  sandy 
bottoms  as  a  function  of  frequency:  (a)  velocity  and  (b)  attenuation. 


25  Hz  by  a  team  from  Lamont-Doherty  Earth  Observatory 
(LDEO)  of  Columbia  University.14  For  comparison,  both 
shear  velocity  and  attenuation  measured  by  the  SAX99 
teams  at  two  frequencies  are  plotted  in  Fig.  2.  Although  only 
one  value  of  shear  velocity  at  1000  Hz  or  25  Hz  from  the 
SAX99  measurements  is  used  as  a  constraining  condition  in 
this  work,  both  the  LF  field-inverted  shear  velocity  and  the 
attenuation  are  close  to  the  SAX99  measurements  at  two  fre¬ 
quencies.  Based  on  Buckingham,15  “the  sole  reported  data 
taken  at  a  frequency  of  1000  Hz”  by  Richardson  from  the 
North  Sea  site  Cl  are  also  plotted  in  Figs.  2(a)  and  2(b). 

For  case  A  (the  velocity  from  the  SAX99  measurements 
at  1000  Hz  as  a  constraining  condition),  the  shear  wave  ve¬ 
locity  and  attenuation,  shown  in  Fig.  2  by  solid  lines,  corre¬ 
spond  to  the  LF  field-inverted  complex  frame  shear  modulus 

H  =  2.91(1  -  0.0981/)  x  107  Pa.  (13) 

In  case  B  (the  velocity  from  the  SAX99  measurements  at 
25  Hz  as  a  constraining  condition),  cs(f)  and  is(f),  shown  in 
Fig.  2  by  dashed  lines,  correspond  to  the  complex  frame 
shear  modulus 
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/<  =  4.46(1  -  0.0426/)  x  107  Pa. 


Seabed  reflection  loss=  8.686Q0  (dB/rad) 


(14) 

Comparing  this  with  Eq.  (13),  a  higher  shear  wave  velocity 
(frame  shear  modulus)  requires  smaller  shear  wave  attenua¬ 
tion  (imaginary  frame  shear  modulus).  This  is  the  shear 
wave  velocity-attenuation  coupling  discussed  before  on 
Eq.  (8). 

By  using  any  set  of  the  shear  wave  velocity  and  attenua¬ 
tion  shown  in  Fig.  2  as  inputs  for  Eq.  (8),  the  LF  field- 
inverted  equivalent  sound  attenuations  a 2e(J)  are  corrected 
to  same  compressional  wave  attenuations,  which  can  be 
expressed  by 

a2  =  (0.420  ±  0.035)  (/'/ 1000) (2  000±0  045)  dB/m 

for  50 -500  Hz.  (15) 

The  original  LF  field  inverted  equivalent  sound  attenuation 
from  Ref.  1,  a2£,  and  shear  wave  corrected  compressional 
wave  attenuation,  a2,  as  well  as  their  power-law  fittings, 
expressed  by  Eqs.  (1)  and  (15),  are  plotted  in  Fig.  3. 

As  an  example  for  comparison,  using  Eqs.  (10) — (12)  we 
calculate  the  near-grazing  bottom  reflection  loss/phase  shift 
parameters  and  the  hidden  depth  of  shallow  water  ( Q ,  P,  and 
AH)  for  two  cases:  with  or  without  shear  waves  that  are  plot¬ 
ted  in  Figs.  4(a)-4(c),  respectively.  With  the  exception  of  the 
shear  modulus  expressed  by  Eq.  (13),  which  is  inverted  in 
this  paper,  other  Biot  parameters  are  the  same  as  those  listed 
in  Ref.  1  (Table  X,  the  column  labeled  “LF  fit”)-  For  the 
case  of  “no  shear  wave,”  we  assume  both  frame  bulk  and 
shear  moduli  are  real  (have  no  imaginary  term).  Figure  4(a) 
shows  strong  frequency  dependence  for  the  near-grazing 
seabed  reflection  loss  from  the  sandy  ocean  bottoms  in  the 
cases  with  or  without  shear  waves.  This  phenomenon  may 
be  extrapolated  from  papers  published  by  Stoll  and  Kan,8 
and  by  Isakson  and  Neilsen.31  An  interesting  frequency  de¬ 
pendence  shown  in  Fig.  4(a)  predicts  that  the  water-sediment 
interface  for  sandy  bottoms  may  act  as  a  filter  with  respect  to 


Average  LF  sound  attenuation  from  1-20  sites 


FIG.  3.  (Color  online)  Comparison  of  the  LF  field  inverted  effective  bottom 
sound  attenuation  ot2E  with  shear  wave  corrected  sound  attenuation  a 2- 


Seabed  reflection  phase  shift  p=n-  Pe 


(b)  Freqency  (Hz) 


FIG.  4.  (Color  online)  Comparisons  between  with  and  without  shear 
waves,  (a)  The  near-grazing  bottom  reflection  loss  factor  Q,  bottom  reflec¬ 
tion  loss  defined  by  201og|V(0)|  =  8.686(20.  (b)  The  near-grazing  bottom 
reflection  phase  shift  parameter  P.  (c)  The  hidden  depth  of  shallow  water 
AH. 
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broadband  long-range  sound  propagation  where  lower-order 
modes  with  low-grazing  angles  are  predominant. 

III.  SUMMARY  AND  DISCUSSION 

The  LF  field  inverted  sound  speed  and  attenuation  in 
sandy  seabottoms  are  “effective”  results  and  yield  a  frequency 
exponent  of  1.8  for  attenuation.1  Whereas  several  physical  the¬ 
ories  on  seabed  acoustics  yield  a  result  of  2.0. 17-20  One  of  pos¬ 
sible  sources  for  this  discrepancy  may  be  neglected  shear 
effects.21  The  LF  field  inversion  results  for  the  exponent  of  1.8 
make  the  assumption  of  an  all-fluid  medium  with  no  shear 
properties.  Based  on  some  theoretical  results  of  Tindle,  Zhang, 
and  Chapman,25--0  the  LF  field-inverted  compressional  wave 
velocity  and  effective  attenuation  of  Ref.  1  are  used  to  estimate 
the  average  shear  wave  velocity  and  attenuation,  and  complex 
shear  modulus  for  a  frequency  band  of  50-500  FIz  in  upper 
layers  of  sandy  and  sand-silt  mixture  ocean  bottoms.  The  re¬ 
sultant  shear  wave  velocity  and  attenuation  are  consistent  with 
the  SAX99  measurement  at  25  Hz  and  1000  Hz. 

Figure  4(c)  shows  that  from  50  to  500  Hz,  the  hidden 
depths  in  shallow  water,  AH,  is  equal  to  about  28.0  m  to 
2.8  m.  The  effective  penetration  depth  of  sound  waves  at 
long  ranges  is  less  than  AH.  Thus,  resultant  shear  wave  ve¬ 
locity  and  attenuation  shown  in  Figs.  2(a)  and  2(b),  as  well 
as  the  complex  shear  modulus  expressed  by  Eqs.  (13)  and 
(14)  can  be  reasonably  considered  as  average  values  in  an 
upper  layer  of  about  20  m  thick. 

In  general,  “shear  wave  properties  are  difficult  to  mea¬ 
sure:  their  attenuation  is  high,  and  it  is  difficult  to  generate 
waves  with  predominantly  transverse  particle  motions.”32 
Buckingham  indicated  “In  situ  measurements  of  shear  wave 
properties  as  functions  of  frequency  are  rarer  than  those  of 
compressional  waves.  In  fact,  no  published  in  situ  measure¬ 
ments  of  shear-wave  dispersion  and  attenuation  over 
extended  frequency  ranges  are  known  to  the  author.”16 
Much  of  the  early  data  on  shear  wave  in  the  literature  from 
both  field  and  laboratory  work  is  in  the  high-frequency 
range.  Most  of  recent  in  situ  measurements  use  interface 
waves  on  the  bottoms;  resultant  shear  wave  speed  and 
attenuation  exhibit  strong  depth  dependence  beneath  the 
water-seabed  interface.  ’  Limited  to  the  authors’ 
knowledge,  very  few  of  those  data  are  in  a  frequency  range 
of  50-1000  Hz  that  we  are  interested  in.  Thus,  the  LF  field 
inverted  shear  wave  velocity  and  attenuation,  and  their  possi¬ 
ble  range  presented  in  this  paper,  might  be  helpful  for  ana¬ 
lyzing  the  effects  of  the  shear  wave  on  LF  long-range  sound 
propagation  in  shallow  water. 

With  the  exception  of  common  uncertainties  that  most 
geoacoustic  inversions  meet,  the  reliability  of  the  LF  field- 
inverted  shear  wave  velocity  and  attenuation  in  this  paper 
mainly  depends  on  the  quality  of  the  LF  field-inverted  effec¬ 
tive  sound  attenuation  in  Ref.  1,  and  the  quality  of  shear 
wave  velocity  measured  by  the  SAX99  group  at  25  Hz  or 
1000  Hz.  It  also  depends  on  the  validity  of  those  seabed  geo¬ 
acoustic  models  that  predict  that  the  compressional  attenua¬ 
tion  varies  as/-2  at  low  frequencies.17-21  Uncertainties  on 
the  LF  field-inverted  effective  sound  attenuation  were 
discussed  in  Ref.  1.  Two  lines  in  Figs.  2(a)  and  2(b)  show 
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variations  of  the  LF  held  inverted  shear  wave  speed  and 
attenuation  caused  by  two  different  constraining  values  of 
shear  wave  velocity  from  two  SAX99  teams  (25  Hz  for 
deeper  sediments,  1000  Hz  for  a  thinner  top  layer).  That  is, 
depending  on  the  constraint  that  we  choose,  we  may  end  up 
with  different  values  for  the  outputs.  However,  these  varia¬ 
tions  (uncertainties)  have  no  effect  on  the  physics  problems 
we  discussed,  including  the  result  shown  in  Fig.  3.  Figures 
2(a)  and  2(b)  and  Eqs.  (13)  and  (14)  show  that  a  broadband 
shear  wave  speed,  attenuation  and  complex  frame  shear 
modulus  may  be  estimated  from  the  LF  field-inverted 
effective  bottom  attenuation  plus  one  value  of  shear  speed 
(or  attenuation)  at  one  single  frequency.  It  will  be  desira¬ 
ble  to  simultaneously  invert  the  effective  bottom  attenua¬ 
tion  from  the  broadband  long-range  held  measurements, 
and  to  measure  shear  wave  velocity  and  attenuation  for 
several  frequencies  by  other  methods  at  the  same  site, 
then  to  have  sufficient  comparison  and  validation  of  the 
results. 
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APPENDIX:  BIOT-STOLL  EXPRESSIONS  ON  SHEAR 
WAVE  VELOCITY  AND  ATTENUATION 

For  a  harmonic  plane  shear  wave  of  cxp(ia>t  —  iksx)  trav¬ 
eling  through  a  sediment  medium,  in  the  Biot-Stoll  model  the 
frequency  equation  of  shear  wave  can  be  expressed  by6'7 

pm2  —  / ik 2  PfCO2 

pfto2  mm 2  —  icoFri/K 

=  (pm2  —  pk2)(mm 2  —  imFy/tc)  —  p2mA  =  0, 

(Al) 

where  m  is  the  angular  frequency,  ks  is  wavenumber,  and  p 
is  the  total  mass  density 

p  =  Ppf  +  (l-0)p,.  (A2) 

Pj  is  the  pore  fluid  mass  density,  rj  is  the  pore  fluid  viscos¬ 
ity,  k  is  the  hydraulic  coefficient  of  permeability,  /f  is  the 
porosity,  ps  is  the  density  of  grain,  and  the  frame  shear 
modulus 

p  =  p0(l  -  ids)  (A3) 

is  a  complex  modulus  to  account  for  dissipation  loss  at  grain 
contact,  and  Ss  is  the  shear  log  decrement.  The  parameter 

m  =  Cjf  (A4) 
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accounts  for  the  phase  of  fluid  flow  with  respect  to  the  macro¬ 
scopic  pressure  gradient.  The  structure  constant  c,  or  tortuosity, 
is  equal  to  1  for  uniform  pores  that  are  parallel  to  the  pressure 
gradient,  and  is  equal  to  3  for  randomly  oriented  pores. 

The  factor  F  is  a  viscosity  correction  to  account  for 
frequency-dependent  viscous  losses  of  oscillating  flow  in  the 
sediment  pore.  Assuming  cylindrical  pores,  Biot  derived  an 
expression  for  the  viscosity  correction  factor  F 


F{£) 


1  eT{e) 

4[l-(2  T{E)/iE)Y 


with 


ber'(s)  +  bei' (a ) 
her(s )  +  bei(e) 


(A5) 


(A6) 


where  her  and  bei  are  the  real  and  imaginary  parts  of  the 
Kelvin  function  of  the  first  kind  of  order  zero,  her'  and  bei' 
are  their  derivatives.  The  argument  s  is 


The  parameter  a  is  the  “pore  size  parameter”  that  depends 
on  the  shape  and  size  of  the  sediment  pores;  it  can  be  expressed 
by 


a  = 


(A8) 


From  Eq.  (A2),  we  may  derive  a  complex  wave  number  for 
shear  wave  in  sediments  as 

k]  =  co2-[l  -  Pf/prri],  (A9) 

F 

where 

,  iFii 

m  =  m - .  (A10) 

KOJ 

Thus,  the  shear  wave  velocity  in  the  sediments  is7,8,10  i7 

cs  =  Re[0/p)1/2(l  -  pj/pm')~1/2}.  (All) 


The  shear  wave  attenuation  is 

as  =  8.686 \m[(p/ p)x^2{  \  —  pj / pm')~1^2}  dB/m.  (A  12) 
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Abstract.  Research  on  ocean  reverberation  has  practical  and  scientific  significance. 
Much  progress  has  been  made  in  the  past  three  decades  to  improve  our  understanding  of 
reverberation.  However,  there  remain  important  unanswered  questions  and  a  real  scarcity  of 
high-quality  basic  research  data  sets.  New  progress  on  the  reverberation  modeling  and  the  low- 
frequency  (LF)  seabed  scattering  characterization  in  shallow  water  (SW)  requires  three  essential 
conditions:  1).  A  reliable  reverberation  model  using  a  physics-based  seabed  scattering  function, 

2).  Carefully  calibrated  broadband  reverberation  data,  and  3).  A  ground  truth  about  the  seabed 
geoacoustic  model.  Some  related  work  on  these  topics  is  introduced  in  this  paper.  The  energy 
flux  method  for  SW  reverberation  is  briefly  introduced.  Integration  of  this  method  with  physics- 
based  seabed  scattering  models  directly  and  intuitively  results  in  a  general  expression  for  SW 
reverberation.  A  simple  relationship  between  the  classic  scattering  cross-section  and  the  modal 
scattering  matrix  is  derived.  Some  basic  research  data  sets,  including  the  reverberation 
level/vertical  coherence  as  a  function  of  time,  frequency,  depth/hydrophone  separation  and  sea 
state,  are  reported.  Reverberation  data  and  model  predictions  are  in  good  agreement,  which 
results  in  some  inversion  results.  The  HF  seabed  scattering  models  and  the  energy  flux  method- 
derived  reverberation  model  are  validated  using  LF  reverberation  broadband  data. 

Keywords:  Reverberation,  seabed  scattering,  energy  flux  method,  modal  scattering  matrix. 

PACS:  43.30.Gv,  43.30.Pc,  43.30.Dr,  and  43.30.Hw 

I.  INTRODUCTION 

Ocean  reverberation  is  often  the  most  severe  limiting  factor  in  relation  to  the  use  of 
active  sonar  systems  in  shallow  water.  Much  progress  has  been  made  to  improve  our 
understanding  of  reverberation,  including  reverberation  modeling  '  and  high 
frequency  (HF)  seabed  scattering. “  “However,  there  remain  important  unanswered 
questions  and  a  real  scarcity  of  high-quality  basic  research  data  sets.”26  SW 
reverberation  involves  seabed  (surface)  scattering  as  well  as  two-way  sound 
propagations,  which  are  controlled  significantly  by  a  geo-acoustic  model  of  bottom. 
New  progress  on  reverberation  is,  in  the  final  analysis,  conditioned  by  our  physical 
understanding  of  both  seabed  (and  surface)  geo-acoustic  models  and  scattering 
models.  Some  related  work  on  these  topics  is  introduced  in  this  paper,  including 
theoretical  modeling,  seabed  geaocoustic/scattering  models,  reverberation  data-model 
comparisons  and  inversions. 


II.  THE  THORETICAL  MODELING:  THE  ENERGY  FLUX 
METHOD  FOR  SW  REVERBERATION 


1  2 

SW  reverberation  modeling  includes  the  virtual  image  method  ,  ray  method  ,  the 
normal-mode  method3-5,  the  energy  flux  method6"15,  the  Green’s  function  and  the 
boundary  perturbation  method  "  ,  the  two-way  coupled  mode  method-  "  ,  the  finite 
element  method23,  the  parabolic  equation  method24,  etc.  This  paper  mainly  introduces 
the  energy  flux  method  for  SW  reverberation. 


A.  General  Reverberation  Expression 


The  energy  flux  (angular  spectrum)  method  for  modeling  SW  reverberation,  based 
on  the  WKB  approximation  of  the  normal-mode  theory,  was  first  presented  in 
archived  Chinese  journals.6"9  Because  it  is  simple  and  intuitive,  the  energy  flux 
method  has  further  been  developed  and  widely  used  for  SW  reverberation 
calculations.10"15 

The  sound  field  in  shallow  water,  generated  by  a  point  source  (with  a  source 
level  =0  dB),  can  be  expressed  as  a  summation  of  nonnal  modes: 

r  2  ni^1  ■ 
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The  sound  intensity  is  expressed  by 
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where  is  the  eigenfunction,  kn  is  the  longitudinal  wave  number,  pn  is  the  modal 
attenuation  factor  of  the  nth  mode,  z0  is  the  source  depth  and  z  is  the  receiver  depth. 


Details  of  the  interference  fine  structures  are  sometimes  of  no  use  for  sonar 
performance  predictions.  Thus,  the  2nd  tenn  of  Eq.  (2)  can  be  neglected. 

According  to  Zhou  6:  (a)  Using  the  W.K.B  approximation  to  the  mode  solution  of 
the  wave  equation;  (b)  Expressing  (slowly  varying)  energy  depth  distribution  of  the 
nth  mode  by 


I  =■ 


(3); 


Sn  tan  0H(z) 

(c)  Decomposing  each  normal-mode  into  a  pair  of  up-  and  down-going  local  quasi¬ 
plane  waves  (Their  grazing  angles  satisfy  kn  =  k(z)  cos  6n  (z)  );  and 


(d)  Using  dn  =  5’;iA:(z0)sin6,H(z0)£/6,„(z0)/2^  to  change  the  incoherent  summation 


over  modes  into  an  integration  over  the  mode-ray  angle  On  ( z0 ) , 

we  can  derive  a  general  expression  for  the  average  field  intensity  in  shallow  water 
given  by: 


(4) 
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Here,  Sn  is  the  cycle  distance  of  the  nth  mode-ray  which  is  related  to  modal  wave 
numbers  by 

Sn=2n/{kn-kn+l)  (5) 

aw  is  the  water  column  absorption  coefficient;  6( z0 )  is  the  grazing  angle  of  the  mode- 
rays  at  the  source  depth  and  6{z)  is  the  grazing  angle  at  the  receiver  depth  (not  just  at 
the  source/receiver  location).  The  modal  attenuation  factor,/?,,  can  be  expressed 
approximately  in  terms  of  the  bottom  reflection  coefficient,  V{6 ) : 

In |  V(#n)  | 


A=- 


(6) 


where  A„  is  a  beam  displacement  on  the  bottom  reflection.  Using  different  methods, 

both  Brekhovskikh"  ’  ,  Smith  ,  Brekhovskikh  and  Lysanov  ,  and  Katsnelson  and 

Petnikov  obtained  an  expression  that  is  similar  to  (4). 

For  calculating  the  average  reverberation  intensity  in  shallow  water,  as 
weighting  functions,  the  two-way  propagation  angular  spectra  and  the  classic  seabed 
scattering  function  (cross  section)  are  well  connected  in  the  angular  domain: 
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Here  h  is  the  water  depth,  <j[6m{h),6n{h)\  can  be  any  seabed  scattering  cross  section 

with  different  scattering  mechanisms  in  consideration.  A  is  a  seabed  scattering  area. 

The  angular  spectrum  of  reverberation  is  expressed  by: 

r21n|F[6> „(/i)]|  n  r21n| V\6jM^ 

exp[ - 1 - -r]  exp[ - 1 - r] 
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B.  A  Reverberation  Expression  in  the  Angular  Domain  for  the 
Pekeris  Waveguide 

For  a  special  case,  in  the  Pekeris  waveguide  (an  iso-velocity  water),  Eq.  (7) 
becomes: 


^  m^rtone. 

R(r)  =  ^h ie-2a"rffe  h  y'Vh(6m,0n)y>e  h  dOmdOn  (9) 
h  r 

Here  the  seabed  scattering  area  A « /ZTcwr0 .  That  is,  the  waveguide  dispersion 
spreading5  has  been  assumed  to  be  much  less  than  the  signal  pulse  duration,  r0 . 


C.  A  General  Reverberation  Expression  in  the  Modal  Domain  for  a 
non-Pekeris  Waveguide 

For  a  non-Pekeris  SW  waveguide,  Eq.  (7)  is  changed  back  to  a  summation  over  the 
mode  numbers,  resulting  in  a  general  expression  in  the  modal  domain: 


R{r,z-z0)  =  ^e-2<v XZ ~ >  °”(z)  £  x  A  x  a{6m,0n)  x  |®m  Wf  |O„(/0f  (10) 


Using  (5)  and  the  property  of  the  slowly  depth-varying  energy  dependence  of  the 


nonnal  mode  |  ®„(z)  |2 ,  Eq.  (10)  can  be  expressed  by: 
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where  D(z)  is  a  small  modification  factor  introduced  by  Zhang  and  Jin  to  overcome  a 
problem  when  a  receiver  is  near  a  turning  point, 


D(z)  =  0.875 


1  dc(z ) 
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(12) 
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Eq.  (11)  is  a  simple  and  practical  expression  for  SW  reverberation  in  the  modal 
domain.  Using  the  classic  cross  sections  of  seabottom  scattering  <j(6m,6n)  and  any 

available  normal-mode  code,  Eq.  (11)  can  be  used  to  calculate  reverberation  intensity 
in  a  SW  waveguide  with  an  arbitrary  velocity  profile.  It  only  requires  two  parameters: 
real  and  imaginary  parts  of  the  modal  eigenvalue,  km  and  fim  . 


D.  Reverberation  Vertical  Coherence  in  Shallow  Water 


Spatial  coherence  in  shallow  water  may  be  defined  as  a  measure  of  arrival  angles’ 
spreading  of  signals.  The  vertical  coherence  of  the  SW  sound  field  and  the  angular 
power  spectrum  of  the  arrival  direction  have  a  relationship  that  is  similar  to  the 
Fourier  transform:7' 32  34 


pv(Az,r,z;z0) 


J  Igps  (&■•  U  z;  zo)e~jkAzsin(0)d  0 
\Iaps(80'’z’zo)d8 


(13) 


Where  A z  is  the  vertical  separation  of  a  pair  of  hydrophones,  and  k  is  the  wave 
number.  Using  a  seabed  scattering  cross-section  and  the  reverberation  angular 


spectrum R  (r,z,z0)  in  (8),  an  expression  for  a  normalized  vertical  cross-correlation 
can  be  obtained  either  in  the  angular  domain  or  in  the  modal  domain.34 


E.  The  Effective  Geoacoustic  Model  for  Sandy  Seabottoms 


Without  a  reliable  geoacoustic  model,  all  data-model  comparisons  in  forward 
physics  problems  of  reverberation,  or  reverberation  inversion  results,  are  questionable. 

The  debate  on  the  sound  speed  dispersion  and  the  frequency  dependence  of  sound 
attenuation  in  seabottoms  has  persisted  for  decades.  One  of  our  recent  papers  analyzes 
and  summarizes  low-frequency  (LF)  measurements  in  shallow  water  that  have  resulted 
in  the  identification  of  nonlinear  frequency  dependence  of  sound  attenuation  in  the 
effective  media  of  sandy  seabottoms  (shown  in  Fig.  1).  Both  the  LF  field-derived 
sound  speed  and  the  attenuation  can  be  well  described  by  the  Biot-Stoll  model,  using 
parameters  that  are  consistent  with  either  theoretical  considerations  or  experimental 
measurements.  ~  Based  on  this  effective  Biot  model  from  the  LF  field  measurements, 
the  bottom  reflection  loss  at  low  grazing  angles  is  shown  in  Fig.  2.  The  interesting 
frequency  dependence  in  Fig.  2  predicts  that  the  water-sediment  interface  for  sandy 
bottoms  may  act  as  a  filter  with  respect  to  the  broadband  long-range  sound  field  where 
lower-order  modes  with  low-grazing  angles  are  predominant.  It  also  affects  the 
frequency  dependence  of  both  long-range  reverberation  and  its  inverted  geophysical 
parameters. 


FIG.  1.  Effective  sound  attenuation  in  sandy 
bottoms  from  LF  field  measurements  0 


Effective  Biot  Model  from  LF  field  measurements 
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FIG.  2.  Bottom  reflection  loss  vs. 

frequency  at  low  grazing  angles 


F.  Seabed  Scattering  Models 

1.  Semi-empirical  Seabed  Scattering  Cross-sections 

For  convenience,  semi-empirical  seabed  scattering  cross-sections  are  often  assumed 
separable  for  incident  angle  and  scattering  angle6'8 

<j(6i,  9S)  =  //sin"/2  6i  sin"72  6s  (14) 

A  special  case  for  n  =  2  is  the  Lambert  law,  which  is  widely  used  for  HF  data- 
model  comparisons.  However,  this  model  has  no  physical  basis. 


2.  Physics-based  Seabed  Scattering  Cross-sections 

In  the  past  30  years,  one  of  the  major  accomplishments  in  ocean  acoustics  is  the 
improvement  in  our  understanding  of  seabed  scattering,  resulting  from  a  significant 
effort,  both  in  at-sea  measurement  and  in  theoretical  modeling,  made  by  many 
investigators.  These  accomplishments,  including  the  rough  bottom  scattering  (RBS) 
model  and  the  sediment  volume  scattering  (SVS)  model,  have  been  summarized  in  an 
underwater  acoustics  monograph  by  Jackson  and  Richardson.  "  According  its 
introduction,  “This  monograph  will  emphasize  small  spatial  scale  and  ‘high 
frequencies’,  very  roughly,  frequencies  from  10  kHz  to  1  MHz”.  We  will  use  SW 
reverberation  data  to  test  their  suitability  at  low  frequencies. 


G.  Modal  Scattering  Matrix  in  SW  Waveguide 


Using  the  energy  flux  method  for  SW  reverberation,  we  have  proved  that,  in  Eq. 
(10)  the  modal  scattering  matrix  in  the  SW  waveguide  between  mode  m  and  mode  n 
can  be  expressed  by: 


MSM(m,n)  =  cr(0m,0n)x\®m(h)\  \On(h)\ 
O2  (h) 

mK)  xfr(0ffl,0„)x 


(15) 


i +U«.)' 


|i  +  r_0U| 

Here  0m  =  cos  '[A',,,  /  kw(h)\ ,  and  (bm(h)  is  the  amplitude  of  the  mth  mode  at  the  water- 
bottom  interface.  V  (0)  is  the  bottom  reflection  coefficient  of  mode-decomposed 


plane  waves.  <j(0m,0n )  can  be  any  physics-based  seabed  scattering  cross  section.  For 
example,  inputting  the  roughness  bottom  scattering  (RBS)  model  into  (15)  can  directly 
result  in  an  expression  for  modal  scattering  matrix  in  the  SW  waveguide  that  is  the 
same  as  that  derived  from  the  Green’s  function  and  the  boundary  perturbation  method 
by  Tracey  and  Schmidt.16,17 


H.  A  Simple  Formula  to  Derive  LF  Seabed  Scattering  at  Low 
Grazing  Angles  in  the  Pekeris  Waveguide 

Low-frequency  seabed  scattering  at  low  grazing  angles  (LGA)  is  almost 
impossible  to  directly  measure  in  shallow  water  (SW),  except  through  inversion  from 
reverberation.  The  closed-form  expressions  for  SW  reverberation  in  the  isovelocity 
waveguide  derived  from  the  energy  flux  method  show  that,  at  a  given  distance  in  the 
three-half  law  region  the  complex,  seabed  scattering  by  all  modes/rays  with  different 
incident  and  scattering  angles  can  equivalently  be  represented  by  the  seabed  back- 
scattering  at  a  single  angle  of  0(r)  =  *Jh  /  (nQr)  .  This  equivalent  angle  depends  on 
distance r ,  seabed  reflection  loss  factor  Q  and  water  depth/? .  The  simple  equivalent 
relationship  can  be  expressed  by: 

BBS(0eg)  =  RL(r)  -10\ogl0(cT0?r2  /  4QHr2)  +  attwx2r  (16)  or 

BBS(0eg )  =  RL(r)  +  2TL(r)  +  6.C )dB  -  101og10  A(r)  +  attw  x  2 r  (17) 


Hereafttvis  the  sound  absorption  in  the  water  column  in  dB/m.  r0is  the  signal 
duration.  A  =  7rrcr0  is  the  seabed  scattering  area.  Eq.  (16)  or  (17)  can  be  used  to 

derive  the  LF&LGA  seabed  back  scattering  strength  by  measuring  the  reverberation 
level  (RL)  or/and  transmission  loss  (TL).  The  method  is  valuable  because  it  allows 
the  seabed  backscattering  strength  to  be  inverted  at  frequencies  and  angles  where  no 
practical  method  of  measurements  exists. 


III.  REVERBERATION  MEASUREMENTS 


A  quality  database  of  reverberation  measurements  is  an  absolutely  essential 
component  of  any  effort  to  understand  the  shallow-water  reverberation  problem. 
However,  to  get  wideband  reverberation  data  for  both  short-  and  long-range  data  in  a 
frequency  range  of  100-3000  Hz  is  a  delicate  task  that  can  be  subject  to  errors. 
Through  years’  efforts,  a  quality  database  of  broadband  reverberation  has  been  set 
up/accumulated.  This  database  has  been  (and  will  be)  used  to  compare  with 
theoretical  reverberation  models  and  to  characterize  the  LF  seabed  scattering.  Some 
examples  selected  from  this  database  are  shown  in  this  section. 

A.  Reverberation  Level  RL(t,f)  as  a  Function  of  Time  and  Frequency 


Reverberation  measurements  were  carefully  calibrated.  RL  was  normalized  to 
the  source  level  (SL).  Fig. 3  shows  RL(t)  of  700-2500  Hz,  obtained  on  June  3,  2012 
at  the  ASIAEX  site  in  the  East  China  Sea.  Fig. 4  shows  RL(f )  at  4  seconds  after  an 
explosive  source  is  detonated,  obtained  from  four  sites. 


Time  ( sec ) 

FIG.  3.  RL(t)  at  the  ASIAEX  site. 
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FIG.  4.  RL{t)  at  4s  from  4  sites 


B.  RL  Depth  Dependence  and  Interference  Phenomenon 

The  first  China-US  joint  ocean  acoustics  experiment  (YS96)  was  conducted  in 
the  summer  of  1996  in  the  Yellow  sea  with  a  strong  thennocline.  The  reverberation 
data  exhibited  some  special  characteristics,  including  strong  depth  dependence,  a 
regular  interference,  a  resonant  volume-scattering  etc.  '  Fig.  5a  shows  the 
RL(t,r;z0)  at  800  Hz  when  a  hydrophone  is  located  below  the  thermocline 


( RD  =  60  m ),  and  two  sources  are  located  below  and  above  the  thennocline, 
respectively  ( SD  =  50 m,  or  7  m).  Fig.  5b  shows  that  the  RL  exhibits  a  regular 
interference  structure  with  increasing  time  (distance),  when  a  receiver  is  in  the 
thennocline  ( RD=22  m)  and  the  source  is  below  the  thermocline  ( SD=50  ). 


Time  (sec)  Time  (  sec  ) 

(a)  (b) 

Fig.  5.  Reverberation  data  at  the  YS96  site  show  (a)  Strong  depth  dependence, 

(b)  Regular  interference  oscillations  when  a  receiver  is  in  the  thermocline. 

C.  Unknown  Resonant  Volume  Scattering  in  Water  Column 

Reverberation  data  at  the  YS96  site  also  showed  an  interesting  resonant 
volume-scattering  in  the  water  column  as  shown  in  Fig.  6:  1).  Around  20-23  seconds 
after  the  explosive  sources  were  denoted  (about  15-16  km  away  from  a  source),  all  5 
shot-produced  reverberation  data  exhibited  strong  water  volume  scattering  between 
700  Hz  and  3000  Hz,  peaked  around  1500  Hz.  It  was  about  7-10  dB  higher  than 
bottom  scattering-induced  RL  from  same  distance.  The  horizontal  size  (dimension)  of 
a  scattering  volume  was  about  2000  meters.  The  reverberation  at  29  km  away  from  the 
source,  produced  by  shot  I,  also  had  a  stronger  wideband  volume  scattering.  But  the 
horizontal  size  of  the  scattering  volume  was  much  smaller.  2)  For  shot  III  (2  minutes 
after  the  first  shot  I),  an  apparent  fish  (group)  produced  resonant  scattering  around  300 
Hz  at  around  29  km  away  from  the  source/receiver  location.  The  scattering  strength  at 
300  Hz  was  about  20  dB  higher  than  the  bottom  scattering  from  same  distance. 
However,  the  wideband  volume  scattering  at  same  distance  shown  in  the  shot  I 
disappeared.  Both  the  source  and  the  receiver  were  located  below  the  thermocline. 
The  scattering  mechanism  that  responded  to  the  observed  volume  scattering  is 
unknown. 

D.  Reverberation  Vertical  Coherence 

Reverberation  vertical  coherence  (RVC)  from  three  sites  has  been  analyzed  and 
added  to  our  quality  data  base  for  further  reverberation  modeling  and  for  the  LF 
seabed  scattering  characterizations.  For  example,  Figs.7a-7b  show  the  RVC  as 
function  of  time,  frequency  and  hydrophone  separation,  obtained  from  the  YS’96  site 
and  the  ASIAEX  site,  respectively. 
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SU=bO(m)  nu=b2(m) 


20  30  40 

Time  (sec) 


20  30  40 

Time  (sec) 

FIG.  6.  Reverberation  data  at  the  YS96  site  show  resonant  scattering  in  the  water 
column  around  20s  to  23  s,  peaked  around  at  1500  Hz 
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FIG.  7.  Reverberation  vertical  coherence  measured  at 

(a)  The  YS96  site  ( Az  =  2.0  m  ),  (b)  At  the  ASIAEX  site  ( Az  =  4.0  m  ). 

E.  Effects  of  Sea  State  on  RL(t,z;z0 )  and  RVC(Az,t,f) 

Wideband  reverberation  measurements  were  repeatedly  conducted  at  a  fixed 
location  in  the  center  of  the  ASIAEX  ECS  site.  Data  were  collected  on  June  3rd  and 
June  5th,  2001  using  the  same  measurement  system.  Sound-speed  profiles  during  the 
two  measurement  periods  were  similar  to  each  other.  However,  wind  speed  and  RMS 
surface-wave  height  (cr)  changed  from  2.74  m/s  and  0.10  m  on  June  3rd  to  7.45  m/s 
and  0.33  m  on  June  5th.  (See  Fig.  8,  where  the  wind  speed  was  measured  at  14  m 
elevation).  Thus,  these  measurements  offer  an  opportunity  to  evaluate  sea  surface 
effects  on  reverberation  level  (RL),  reverberation  vertical  coherence  (RYC),  and 


RVC-inverted  bottom  acoustic  parameters  in  shallow  water  (SW). 40,41  Two  sets  of  data 
for  RVC  and  RL  in  a  frequency  range  of  100-2500  Hz  show  differences  that  are  the 
apparent  effects  of  the  surface  roughness.  In  contrast  to  many  peoples’  intuition:  with 
increasing  sea  state,  the  RVC  increases  (shown  in  Fig.  9),  but  the  RL  decreases 
(shown  in  Fig.  10). 


surface  data,  APL-UW 


FIG.  8.  Wind  speed  and  surface  wave  height  variability  during  the  ASIAEX  experiment. 
(Collected  and  analyzed  by  Peter  Dahl,  APL-UW) 


F  =  800  Hz,  Az=1.7m,  ASIAEX  site 
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FIG.  9.  The  RVC  increases  with  increased  sea  state,  (a)  800  Hz,  (b)  1200  Hz 


F  =  1250  Hz,  SD  =  50m  ,  RD  =  56.5m 


F  =  2500  Hz,  SD  =  50  m  ,  RD  =  56.5m 


FIG.  10.  The  RL  decreases  with  increased  sea  state,  (a)  1250  Hz,  (b)  2500  Hz 


IV.  REVERBERATION  DATA-MODEL  COMPARISON  AND 

INVERSION 


A.  Data-model  Comparisons  Using  Physics-based  Seabed  Scattering 
Model 


Two  reverberation  measurements  were  conducted  in  the  East  China  Sea.  Two 
sound  velocity  profiles  (SVP)  are  shown  in  Fig.  1 1 .  The  SVP  at  the  ECS  site  I  was 
almost  iso-velocity  (the  Pekeris  shallow  water).  The  SVP  at  the  ASIAEX  site  on  June 
3,  2001  had  a  weaker  negative  gradient. 
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FIG.  11.  Sound  velocity  profiles  at  two  sites  during  reverberation  measurements 
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The  effective  Biot  seabed  geoacoustic  model  derived  from  the  LF  field 
measurements  is  used  in  our  reverberation  modeling.  The  physics-based  rough  seabed 
scattering  cross  section  and  the  sediment  volume  cross  section  have  been  extensively 
tested  with  FIF  data  ,  and  are  used  in  Eq.  (9)  and  Eq.  (11)  for  data-model 
comparisons.  Figs.  12a  and  12b  show  the  reverberation  data-model  comparisons  at 
the  ECS  site  I,  using  Eq.  (9)  in  the  angular  domain  for  a  Pekeris  waveguide.  Figs.  13a 
and  13b  show  the  data-model  comparisons  at  the  ASIAEX  site,  using  Eq.  (11)  in  the 
modal  domain  for  the  non-Pekeris  waveguide.  All  required  spectral  strength  co2  and 

the  exponent  y2  for  the  bottom  roughness  spectrum,  and  the  sediment  volume 
scattering  cross-section  <jv  for  the  sediment  volume  scattering  are  shown  in  these 

figures.  The  data  and  the  model  predictions  are  in  good  agreement.  That  is,  the 
physics-based  HF  seabed  scattering  models  can  be  used  to  predict  LF  reverberation 
using  a  right  set  of  seabed  geoacoustic  and  scattering  parameters.  Our  reverberation 
expressions  in  both  the  angular  and  the  modal  domain  are  also  validated. 
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FIG.  12.  Data-model  comparison  at  the  ECS  site  I  (the  Pekeris  waveguide)  using  Eq.  (9), 

(a)  315  Hz,  (b)  800  Hz 
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FIG.  13.  Data-model  comparison  at  the  ASIAEX  site  (non-Pekeris  waveguide),  using  Eq.  (11) 

(a)  800  Hz,  (b)  1500  Hz 


B.  Low  Grazing  Angle  Scattering  at  Low  Frequencies36 


The  equivalent  relationship  between  the  multi-mode/ray  scattering  and  one 
specific  angle  backscattering  (Eq.  16)  and  the  LF  field  derived  effective  Biot  seabed 
geoacoustic  model  are  used  to  extract  the  bottom  backscattering  strength  as  a  function 
of  angle  and  frequency.  Figs.  14a  and  14b  show  the  RL  inverted  bottom  backscattering 
strength  as  a  function  of  grazing  angle  and  frequency.  Apparently,  the  FF  bottom 
scattering  exhibits  strong  frequency  dependence.  At  the  small  grazing  angles  of 
4  ,6  ,8°  andl0° ,  the  bottom  backscattering  strength  between  200  Hz  and  2500  Hz  can 
approximately  be  expressed  by  BBS(B)  =  -34.7  +  101og10(/  / 1000)27  ( dB ) . 


For  grazing  angles  of  4,6,8,  Ilf 


Grazing  angle  (  degree ) 

(a) 


(b) 


FIG.  14.  RL  inverted  LF  bottom  backscattering  strength  vs.  (a)  angle,  (b)  frequency 

C.  RVC -inverted  bottom  acoustic  parameters34'41 


Reverberation  vertical  coherence  has  been  used  to  derive  effective  seabottom 
sound  speed,  attenuation  and  reflection  loss.  Fig.  15a  shows  the  best  match  between 
the  measured  and  theoretical  RVC  at  400  Hz  for  the  ASIAEX  site  using  a  pair  of 
values  for  bottom  sound  velocity  and  attenuation.  Fig.  15b  shows  the  effective 
boundary  (bottom  and  surface)  reflection  loss  difference  at  the  ASIAEX  site,  inverted 
from  RVC  on  3  June  and  5  June,  respectively.  It  shows  that  higher  boundary  loss 
corresponds  to  the  increase  of  the  sea  state.  This  figure  also  shows  that  the  additional 
surface  loss  in  shallow  water  is  less  than  the  coherent  mean  reflection  loss  predicted 
by  the  perturbation  theory  with  the  Pierson-Moskowitz  model,  because  some  forward 
scattered  energy  returns  to  the  waveguide.  These  results  were  obtained  using  the  semi- 
empirical  seabed  scattering  model.  In  the  future,  it  is  desirable  to  use  physics-based 
seabed  scattering  models  to  get  better  results,  including  the  sea  surface  forward 
scattering. 
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Boundary  loss  difference  for  two  sea  states 
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FIG.  15.  RVC-inverted  seabottom  parameters  at  the  ASIAEX  site,  (a)  Bottom  velocity  and  attenuation 
at  400  Hz,  (b)  RVC-inverted  boundary  loss  difference  for  two  sea  states  on  3  and  5  June,  2001. 


D.  Reverberation-inverted  LF  seabed  scattering  parameters 

The  reverberation  data-model  comparisons  in  Sec.  IV  (A)  have  proved  that  the 
well  developed  HF  physics-based  seabed  scattering  model  can  directly  be  used  to 
predict  LF  SW  reverberation.  However,  the  broadband  reverberation-inverted  LF 
seabed  scattering  parameters  are  different  from  those  values  often  used  for  HF 
modeling  or  testing.  For  example,  the  HF  roughness  spectrum  exponent  is  restricted 
to  the  range  of  with  a  mean  of  3.0.  But,  the  long-range  reverberation  inverted  LF 
roughness  spectrum  exponent  has  a  much  smaller  value  with  a  mean  of  1 .3.  The  HF 
sediment  volume  scattering  cross  section  ov  is  generally  assumed  to  have  linear 
frequency  dependence.  However,  the  broadband  LF  long-range  reverberation-inverted 
<7V  exhibits  much  stronger  frequency  dependence,  up  to  the  fourth  power  of  the 
frequency. 


V.  SUMMARY 

The  energy  flux  method  for  SW  reverberation  based  on  the  WKB 
approximation  of  normal-mode  theory,  first  presented  three  decades  ago  in  archived 
Chinese  journals,  is  briefly  introduced  with  new  results.  Benefiting  from  the 
accomplishments  of  the  seabed  scattering  research  community,  we  have  integrated  the 
energy  flux  method  for  SW  reverberation  with  the  physics-based  seabed  scattering 
models,  and  obtained  a  general  expression  for  SW  reverberation  in  the  angular  domain 
(Eq.  7)  or  in  the  modal  domain  (Eq.  11).  The  integration  also  directly  results  in  a 
simple  relationship  between  the  classic  boundary  scattering  cross-section  and  the 
modal  scattering  matrix  in  SW  waveguide  expressed  by  Eq.  (15).  The  development  of 
a  high  quality  reverberation  data  base  has  become  a  pressing  task  for  SW 
reverberation  modeling  and  for  LF  seabed  scattering  characterization.  Some  data  on 
RL  and  RVC  as  a  function  of  time,  frequency,  depth/hydrophone  separation  and  sea 
state,  are  reported.  The  reverberation  data  and  the  energy  flux  method-derived 
reverberation  model  are  in  good  agreement.  Reverberation-based  inversion  results, 
briefly  reported  in  this  paper,  include  the  low  grazing  angle  seabed  scattering  strength 
vs.  angle  and  frequency,  seabottom  velocity  and  attenuation,  additional  sea  surface 
reflection  loss,  and  the  LF  seabed  scattering  parameters  for  the  rough  seabed  scattering 
and  the  sediment  volume  scattering  section,  etc. 


ACKNOWLEDGMENTS 


This  work  was  supported  by  ONR  OA322  and  CAS. 


REFERENCES 


1 .  K.  V.  Mackenzie,  “Long-range  shallow -water  bottom  reverberation,”  J.  Acoust.  Soc.  Am.  34,  62-66 
(1962). 

2.  C.  Y.  Wu,  ACTA  ACUSTICA  (Chinese)  4,  114-119  (1979)  C.Y.  Wu,  “Calculation  of  shallow  water 
reverberation  intensity  based  on  ray  theory.  Part  (1),”  in  Acta  Acustica  (Chinese)  4,  114-119  (1979). 
Part  (2)  in  Acta  Acustica  (Chinese),  5,  210-220  (1980). 

3.  H.P.  Bucker  and  H.E.  Morris,  “Normal-mode  reverberation  in  channels  or  Ducts,”  J.  Acoust.  Soc. 
Am.  44,  827-828  (1968). 

4.  R.H.  Zhang  and  G.L.  Jin,  “Normal-mode  theory  of  average  reverberation  intensity  in  shallow-water,” 
J.  Sound  Vib.,  119,  215-223  (1987). 

5.  D.D.  Ellis,  “A  shallow-water  normal-mode  reverberation  model,”  J.  Acoust.  Soc.  Am.  97,  2804- 
2814  (1995). 

6.  J.X.  Zhou,  “The  Analytical  method  of  angular  power  spectrum,  range  and  depth  structure  of  the 
Echo-to-Reverberation  Ratio  in  Shallow  Water  Sound  Field,”ACTA  ACUSTICA  (Chinese),  5  (2), 
86-99  (1980). 

7.  J.X.  Zhou,  “The  vertical  coherence  of  long-range  reverberation  in  the  shallow  water  homogeneous 
layer,”  ACTA  OCEAN OLOGI A  SINICA,  1  (2),  212-218  (Dec.  1979)  (in  Chinese). 

8.  J.  X.  Zhou,  D.H.  Guan,  E.C.  Shang  and  E.S.  Luo,  “Long-range  reverberation  and  bottom  scattering 
strength  in  shallow  water,”  Chin.  J.  Acoust.  1,  54-63  (1982). 

9.  J.X.  Zhou,  X.Z.  Zhang  and  E.S.  Luo,  “Shallow-water  reverberation  and  small  angle  bottom 
scattering,”  in  <Shallow-water  Acoustics>  edited  by  R.H.  Zhang  and  J.X.  Zhou,  pp.  315-322 
(China  Ocean  Press,  1997) 

10.  C.H.  Harrison, “Closed-form  expressions  for  ocean  reverberation  and  signal  excess  with  mode 
stripping  and  Lambert’s  law,”  J.  Acoust.  Soc.  Am.  114,  2744-2756  (2003). 

11.  C.H.  Harrison,  “Closed  form  bistatic  reverberation  and  target  echoes  with  variable  bathymetry  and 

sound  speed,”  IEEE  J.  Ocean  Eng.  30,  660-675  (2005). 

12.  C.W.  Holland,  “Constrained  comparison  of  ocean  waveguide  reverberation  theory  and 
observations,”  J.  Acoust.  Soc.  Am.  120,  1922-1931  (2006). 

13.  M.  A.  Ainslie,  “Observable  parameters  from  multipath  bottom  reverberation  in  shallow  water,”  J. 
Acoust.  Soc.  Am.  121,  3363-3376  (2007). 

14.  C.W.  Holland,  “Fitting  data,  but  poor  predictions:  Reverberation  prediction  uncertainty  when 
seabed  parameters  are  derived  from  reverberation  measurements,”  J.  Acoust.  Soc.  Am.  123,  2553- 
2562  (2008). 

15.  M.  A.  Ainslie,  Principles  of  Sonar  Performance  Modeling,  (Springer-Verlag  Berlin  Heidelberg 
2010),  pp.  495-510. 

16.  B.H.  Tracey  and  H.  Schmidt,  “Seismo-acoustic  field  statistics  in  shallow  water,”  IEEE  J.  Ocean. 
Eng.  22,  317-331  (1997)]. 

17.  H.C.  Song,  S.  Kim,  W.S.  Hodgikiss  and  W.A.  Kuperman,  “Environmentally  adaptive  reverberation 
nulling  using  a  time  reversal  mirror,”  J.  Acoust.  Soc.  Am.,  116  (2),  762-768  (2004). 

18.  E.C.  Shang,  T.F.  Gao  and  J.R.  Wu,  “A  shallow-water  reverberation  model  based  on  perturbation 
theory,”  IEEE  J.  Oceanic  Eng.,  33  (4),  451-461  (2008). 

19.  A.  Galinde,  N.  Donabed,  M.  Andrews,  S.  Lee,  N.C.  Markris  and  P.  Ratial,  “Range-dependent 
waveguide  scattering  model  calibrated  for  bottom  reverberation  in  a  continental  shelf 
environment,”  J.  Acoust.  Soc.  Am.  123,  1270-1281  (2008). 

20.  D.J.  Tang  and  D.R.  Jackson,  “Application  of  small-roughness  perturbation  theory  to  reverberation 
in  range-dependent  waveguide,”  J.  Acoust.  Soc.  Am.  131,  4428-4441  (2012). 

21.  S.  A.  Stotts  and  R.A.  Koch,  “Rough  surface  scattering  in  a  Born  approximation  from  a  two-way 
coupled-mode  formalism,”  J.  Acoust.  Soc.  Am.  125,  EL242-EL248  (2009). 

22.  S.A.  Stotts,  D.P.  Knobles  and  R.A.  Koch,  “Scattering  in  a  Pekeris  waveguide  from  a  rough  bottom 
using  a  two-way  coupled  mode  approach,”  J.  Acoust.  Soc.  Am.  129,  EL172-178  (2011). 

23.  M.J.  Isakson  and  N.P.  Chotiros,  “Finite  element  modeling  of  reverberation  and  transmission  loss  in 
shallow  water  waveguides  with  rough  boundaries,”  J.  Acoust.  Soc.  Am.  129,  1273-1279  (2011). 

24.  J.F.  Lingevich  and  K.D.  LePage,  “Parabolic  equation  simulations  of  reverberation  statistics  from 
non-Gaussian-distributed  bottom  roughness,”  IEEE  J.  Ocean.  35,199-208  (2010). 


25.  D.R.  Jackson  and  M.D.  Richardson,  High-frequency  Seafloor  Acoustics,  (Springer,  New  York, 
2007). 

26.  J.R.  Preston,  “Report  on  the  1999  ONR  Shallow-water  Reverberation  Focus  Workshop,”  Tech. 

Memorandum  No.  99-155,  Applied  Research  Laboratory,  Penn  State  University,  31  Dec.,  1999. 

27.  L.M.  Brekhovskikh,  “The  Average  Field  in  an  Underwater  Sound  Channel,”  Sov.  Phys.  Acoust., 
11(2),  126-134  (1965). 

28.  L.M.  Brekhovskikh,  “Elements  of  sound  field  theory  in  the  ocean,”  in  <Ocean  Acoustics>,e d.  By 
LM.  Brekhovskikh  ,  pp. 79-162,  Nauka,  Moscow  1974  (in  Russian) 

29.  P.M.  Smith,  “Average  sound  transmission  in  range-dependent  channels,”  J.  Acoust.  Soc.  Am.  55, 

1197-1204  (1974). 

30.  L.M.  Brekhovskikh  and  Y.P.  Lysanov,  <Fundaments  of  Ocean  Acoustics>,3rd  Edition,  (Springer- 

Verlag,  New  York,  2003)  Sec.  5.6 

31.  B.  G.  Katsnelson  and  V.G.  Petnikov,  <Shallow-Water  Acoustics>,  (Springer-Praxis,  New  York, 

2002),  Sec.  3.5. 

32.  P.M.  Smith,  “Spatial  coherence  in  multipath  or  multimodal  channels, ”J.  Acoust.  Soc.  Am.  60,  305- 
310  (1976). 

33.  J.X.  Zhou,  “Vertical  coherence  of  the  sound  field  and  boundary  loss  in  shallow  water,”  Chinese 
Phys.,  1,  494-504(1981). [  Acta  Acust.,  4  ,34-43  (1979)] 

34.  J.X.  Zhou,  X.Z.  Zhang,  P.H.  Rogers,  J.A.  Simmen,  P.  H.  Dahl,  G.L.  Jin  and  Z.H.  Peng, 
“Reverberation  vertical  coherence  and  seabottom  geoacoustic  inversion  in  shallow  water,”  IEEE  J. 

Ocean  Eng.  29,  988-999  (2004). 

35.  J.X.  Zhou,  X.Z. Zhang,  and  D.P.  Knobles,  “Low-frequency  geoacoustic  model  for  the  effective 
properties  of  sandy  seabottoms,”  J.  Acoust.  Soc. Am.,  125,  2847-2866  (2009). 

36.  J.X.  Zhou  and  X.Z.  Zhang,  “Low  frequency  seabed  scattering  at  low  grazing  angles,”  J.  Acoust. 
Soc.Am.,  131,  2611-2621  (2012). 

37.  J.X.  Zhou  and  X.Z.  Zhang,  “Shallow-water  reverberation  level:  measurement  techniques  and  initial 
reference  values,”  IEEE  J.  Ocean  Eng.  30,  832-842  (2005). 

38.  J.X.  Zhou,  X.Z.  Zhang,  and  G.L.  Jin,  “Reverberation  characteristics  in  shallow  water  with  a  strong 
thermocline,”  J.  Acoust.  Soc.  Am.  107,  2920  (2000). 

39.  F.  Li,  G.  Jin,  R.  Zhang,  J.  Liu,  and  L.  Guo,  “An  oscilation  phenomenon  of  reverberation  in  shallow 
water  with  a  thermocline,”  J.  Sound  Vib.,  252,  457-468(2002). 

40.  F.H.  Li,  J.J.  Liu  and  R.H.  Zhang,  “A  model/data  comparison  for  shallow-water  reverberation,” 
IEEE  J.  Ocean  Eng.  29,  1060-1066  (2004). 

41.  J.X.  Zhou,  X.Z. Zhang  Z.  H.  Peng  and  J.S.  Martin,“Sea  surface  effects  on  shallow-water 
reverberation,”  J.  Acoust.  Soc.Am.,  121,  98-107  (2007). 


CrossMark 


Integrating  the  energy  flux  method  for  reverberation  with 
physics-based  seabed  scattering  models:  Modeling  and  inversion 


Ji-Xun  Zhoua)  and  Xue-Zhen  Zhang 

School  of  Mechanical  Engineering,  Georgia  Institute  of  Technology ,  Atlanta,  Georgia  30332-0405 
(Received  23  August  2012;  revised  29  March  2013;  accepted  26  April  2013) 

During  the  past  30  years,  one  of  the  major  accomplishments  in  ocean  acoustics  is  the  improvement 
of  understanding  seabed  scattering,  resulting  from  a  significant  effort  of  both  at-sea  measurement 
and  theoretical  modeling.  [Jackson  and  Richardson:  High-Frequency  Seafloor  Acoustics,  2007]. 
Benefiting  from  this  accomplishment,  this  paper  integrates  the  energy  flux  method  for  shallow- 
water  (SW)  reverberation  [Zhou,  (Chinese)  Acta  Acust.  5,  86-99  (1980)]  with  the  physics-based 
seabed  scattering  models.  This  integration  directly  and  intuitively  results  in  general  expressions  for 
SW  reverberation  in  the  angular  and  modal  domains.  The  latter  expression  is  the  same  as  the  modal 
reverberation  expression  derived  from  the  Green’s  function  and  boundary  perturbation  method  by 
Tracey  and  Schmidt  [IEEE  J.  Ocean.  Eng.  22,  317-331(1997)].  The  integration  also  results  in  a 
simple  relationship  between  the  classic  boundary  scattering  cross  sections  and  the  modal  scattering 
matrix  in  SW  waveguides.  The  bottom  roughness  spectrum  and  sediment  volume  scattering  cross 
section  at  low  grazing  angles  are  inverted  in  a  frequency  range  of  150-2500  Hz  from  the  wideband 
long-range  reverberation  data  by  using  the  Biot  seabed  geoacoustic  model.  The  results  may  offer 
some  reference  data  sets  for  future  analysis  of  the  low-frequency  seabed  scattering  mechanisms. 
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I.  INTRODUCTION 

Acoustic  reverberation  in  shallow-water  (SW)  wave¬ 
guides  has  recently  become  a  research  area  of  intensely 
increasing  interest.  Much  progress  has  been  made  in  two 
areas:  SW  reverberation  modeling  and  seabed  scattering. 
However,  these  two  research  communities  have  not  yet  had 
enough  communication  with  each  other. 

SW  reverberation  modeling  includes  the  virtual  image 
method,1  the  ray  method,2'3  the  normal-mode  method,1  8  the 
energy  flux  method,9-20  the  Green’s  function  and  boundary 
perturbation  method,2 1-28  the  two-way  coupled  mode 

no  in  , 

method,  the  finite  element  method,  the  parabolic  equa¬ 
tion  method,32  etc.  However,  most  of  these  models  have  not 
yet  fully  taken  advantage  of  the  progress  made  by  the  seabed 
scattering  community. 

In  general,  seabed  scattering  is  a  dominant  mechanism 
of  generating  SW  reverberation.  In  the  past  30  years,  one  of 
the  major  accomplishments  in  ocean  acoustics  is  the 
improvement  in  our  understanding  of  seabed  scattering, 
resulting  from  a  significant  effort  of  both  at-sea  measure¬ 
ment  and  theoretical  modeling  made  by  many  investigators. 
These  accomplishments,  including  the  rough  bottom  scatter¬ 
ing  (RBS)  model  and  the  sediment  volume  scattering  (SVS) 
model,  have  been  summarized  in  an  underwater  acoustics 
monograph  by  Jackson  and  Richardson.3’  However,  there 
remains  a  real  scarcity  of  seabed  scattering  data  sets  in  the 
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low  frequency  (LF)  range  of  100-2600  Hz  that  can  be  used 
to  compare  with  the  RBS  and  SVS  models. 

Motivated  by  the  above-mentioned  three  “howevers,” 
this  paper  integrates  the  energy  flux  method  of  SW  reverber¬ 
ation  with  the  physics-based  seabed  scattering  models.  The 
integration  directly  and  intuitively  results  in  a  practical 
expression  for  long-range  SW  reverberation  in  terms  of 
seabed  physical/scattering  parameters. 

The  present  paper  first  introduces  a  theoretical  model  for 
SW  reverberation  using  the  energy  flux  method,  then  charac¬ 
terizes  the  LF  seabed  scattering  at  low  grazing  angles  from 
broadband  reverberation  data  using  the  Biot  geoacoustic 
model.  The  paper  is  organized  as  follows.  Section  II  briefly 
introduces  the  energy  flux  (angular  spectrum)  method  for  SW 
reverberation,  which  is  based  on  the  Wentzel-Kramers- 
Brillouin  (WKB)  approximation.9  For  readers’  convenience, 
two  physics-based  seabed  scattering  models  that  are  widely 
used  for  high  frequency  (HF)  seafloor  scattering,  the  RBS 
model  and  the  SVS  model,3’  are  briefly  introduced  in  Sec. 
III.  Section  IV  integrates  the  energy  flux  method  of  SW 
reverberation  with  these  physics-based  RBS  and  SVS  mod¬ 
els.  The  integration  results  in  a  practical  expression  for  SW 
reverberation  in  the  angular  domain  or  in  the  mode  domain. 
A  simple  relationship  between  the  bottom  scattering  cross 
sections  and  the  modal  scattering  matrix  in  the  SW  wave¬ 
guides  obtained  from  the  integration  is  described  in  Sec.  V. 
For  validating  the  energy  flux  method-derived  reverberation 
expressions  with  the  HF  RBS  and  SVS  models,  data-model 
comparisons  for  the  ASIAEX  site  are  given  in  Sec.  VI. 
Section  VII  characterizes  the  LF  seabed  scattering  at  low 
grazing  angles  using  broadband  long-range  reverberation 
data  and  the  Biot  geoacoustic  model.  Finally,  Sec.  VI  sum¬ 
marizes  and  discusses  the  results  of  this  paper. 
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II.  THE  ENERGY  FLUX  METHOD  FOR  SW 
REVERBERATION 


The  angular  spectrum  (energy  flux)  method  for  model¬ 
ing  SW  reverberation,  based  on  the  WKB  approximation  of 
normal-mode  theory,  was  first  presented  in  archived  Chinese 
journals.9-1’  Due  to  its  simplicity  and  intuitiveness,  the 
energy  flux  method  has  further  been  developed  and  widely 
used  for  SW  reverberation  calculations.9-20  This  method  is 
briefly  introduced  here,  and  the  interested  reader  is  referred 
to  the  basic  concepts  in  Sec.  II  and  Sec.  IV  of  Ref.  9. 

The  sound  field  in  SW,  generated  by  a  point  source 
(with  a  source  level,  SL  =  0  dB),  can  be  expressed  as  a  sum¬ 
mation  of  normal  modes34 
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where  ®„  is  the  normalized  eigenfunction,  kn  is  the  longitu¬ 
dinal  wave  number,  (>„  is  the  modal  attenuation  factor  of  the 
uth  mode  (i.e.,  the  complex  wave  number  Kn  =  kn  +  i[ln).  z0 
is  the  source  depth  and  z  is  the  receiver  depth.  The  first  term 
of  Eq.  (2)  is  an  incoherent  summation  of  normal  modes.  It 
describes  the  average  intensity  of  the  sound  field.  For  many 
practical  applications  such  as  sonar  performance  predictions, 
an  average  estimate  of  the  sound  intensity  is  quite  sufficient, 
i.e.,  the  second  term  of  Eq.  (2),  which  describes  the  interfer¬ 
ence  fine  structure,  can  be  neglected. 

According  to  the  WKB  approximation,  the  smoothed 
(slowly  varying)  energy  depth  distribution  of  the  nth  mode 
can  approximately  be  expressed  as 


where  and  rjn  are  the  upper  and  lower  turning  (or  reflec¬ 
tion)  depths  of  the  mh  mode,  and  are  the  correspond¬ 
ing  phase  shifts,  and  kw(z)  =  co/c(z)  is  the  wave  number  in 
water.  The  modal  attenuation  factor,  [!„,  can  be  expressed 
approximately  in  terms  of  the  bottom  reflection  coefficient 
V(0), 
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Each  normal  mode  can  be  decomposed  into  a  pair  of  up- 
and  down-going  local  quasi-plane  waves.  We  shall  call  the 
local  plane  wave  an  equivalent  “mode  ray.”  The  angle  of  the 
mode  ray  and  the  eigenvalue  of  the  /?th  mode  are  related  by 

k„  =  kw(z)  cos  6 „(z).  (7) 


It  is  well  known 
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Using  Eqs.  (7)  and  (8),  we  have 
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The  incoherent  summation  over  modes  can  be  changed  to  an 
integration  over  the  mode-ray  angle  6.  Proceeding  from  Eqs. 
(2),  (3),  (6),  (8),  and  (9),  and  adding  a  term  of  water  absorp¬ 
tion,  we  can  derive  a  general  expression  for  the  average  field 
intensity  in  SW  given  by 
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Here,  S„  is  the  cycle  distance  of  the  mh  mode  that  is  the  sum 
of  the  geometric  cycle  distance  Sng  and  the  beam  displace¬ 
ment  A,,:35 


S 


n 


S„g  +  A„  —  2 


dz 

L  tan  6„(z) 


+  A„, 


(4) 


and  6n  is  the  mode  angle.  A„  is  a  beam  displacement  on 
the  bottom  reflection.  If  the  bottom  reflection  is  expressed 
by  V(8)  =  |V(0)|  exp(i<pfc),  with  a  phase  of  (ph,  then 
A„  =  d(ph/dkn.  The  eigenvalue  k„  of  the  nth  mode  satisfies 
relationship 


K2dz  +  B(n  +  £fjn  =  2nn, 


(5) 


Using  different  methods,  Brekhovskikh, 36,37  Smith,38 
Brekhovskikh  and  Lysanov,  ’9  and  Katsnelson  and  Petnikov40 
obtained  an  expression  that  is  similar  to  Eq.  (10).  In  Eq.  (11), 
ocw  is  the  water  column  absorption  coefficient;  (Hz0)  is  the 
grazing  angle  of  the  mode  rays  at  the  source  depth,  and  ()(z) 
is  the  grazing  angle  at  the  receiver  depth  (note:  Here  “depth” 
means  in  the  whole  water  column  not  just  at  the  source/ 
receiver  location).  The  quantity  defined  by  Eq.  (11),  Iaps,  is 
referred  to  as  the  “angular  power  spectrum”  of  sound  propa¬ 
gation.9  Smith  also  called  it  the  “energy  density  in  direction 
of  arrival”  (Ref.  38).  The  average  sound  intensity  in  SW 
depends  only  on  water  absorption  and  two  other  parameters: 
The  bottom  reflection  coefficient  V(0)  and  the  cycle  distance 
S„  of  the  mode  ray. 

Equation  (10)  is  a  simple  and  intuitive  expression  for 
calculating  averaged  SW  field  characteristics  such  as  sound 
propagation,  reverberation,  noise,  and  their  spatial 
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coherence.  It  is  only  an  angular  weighting  process.  For  cal¬ 
culating  the  average  reverberation  intensity  in  SW,  the  two- 
way  propagation  angular  spectra  and  the  classic  seabed  scat¬ 
tering  function  (cross  section)  are  naturally  connected  in  the 
angular  domain. 

Using  the  single  scatter  approximation,  the  average 
reverberation  intensity  R  can  be  expressed  in  terms  of  Iaps 
by9 


where 


G  = 


[  cos  0 i  cos  8S  cos 


B  = 


sin  ( 


;  Sin  I 


'ps 


a2pap 


(16) 

(17) 
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hr 


=  7e 


-2awr 


exp({2  \n\V[9m(h)}\/Snt}r) 


Sm  tan  8m(h) 

x  Aa[6m(h),9n(h)\ 
exp({21n|V,[0n(A)]|/5„}r) 


S„  tan0„(z) 


dOm(z0)dOn(h ), 

(12) 


where  Iaps{9m:  r,  /;;  zo)  is  the  incident  angular  power  spec¬ 
trum  at  a  distance  of  r,  incident  upon  a  seabed  scattering 
area  A;  Iaps(9n,r,z-,h)  is  the  angular  power  spectrum 
returned  to  a  receiver  at  a  depth  of  z  from  the  same  bottom 
scattering  area;  h  is  a  water  depth;  ab[6m(h),9n(h)\  is  the 
bottom  scattering  cross  section  per  unit  area;  9m(Ji)  is  the 
incident  grazing  angle  of  mode  m  at  the  bottom;  and  9n(h)  is 
the  scattering  angle  of  mode  n. 


III.  TWO  PHYSICS-BASED  SEABED  SCATTERING 
MODELS 

The  bistatic  seabed  scattering  models  for  high  frequen¬ 
cies  have  extensively  been  investigated,  tested  against  data 
with  good  results,  and  summarized  in  a  comprehensive 
monograph  (Ref.  33).  These  HF  scattering  models  have  been 
extended  to  the  geoacoustic  bottom  interaction  model 
(GABIM)  for  LF  range.41  For  the  simplicity  on  the  reverber¬ 
ation  data-model  comparisons,  the  present  study  will  use 
two  (commonly  adopted)  physics-based  models:  The  RBS 
model  and  the  SVS  model  described  in  Ref.  33.  For  readers’ 
convenience,  all  notations  in  the  expressions  for  these  two 
scattering  models  in  present  paper  will  be  same  as  those 
used  by  Jackson  and  Richardson.” 

A.  The  RBS  model 

The  bistatic  scattering  cross  section  for  the  RBS  model 
in  the  small-roughness  perturbation  approximation  is33 


sin  9 pi  =  —  a2  cos2  0,-, 


sin  8ps  =  J 1  -  a2  cos2  6S, 


1  +  iS„ 


pCw  ln(10)  f  cb  x  y.p 


40ti/ 


\2nf  x  8.686/’ 


(18) 

(19) 

(20) 

(21) 


where  W(AK)  is  the  roughness  spectrum,  which  is  the 
Fourier  transform  of  the  interface  relief  covariance. 
Equation  (14)  is  one  of  the  simplest  and  most  often  used  iso¬ 
tropic  forms,  the  “power  law.”  AK  is  the  magnitude  of  the 
two-dimensional  “wave  vector.”  The  parameters  of  the  spec¬ 
trum  are  the  “spectral  strength”  cut  and  the  “spectral 
exponent”  y2.  The  definition  of  angular  coordinates  used  in 
these  equations  is  shown  in  Fig.  2.4  of  Ref.  33,  0,  and  0,  are 
the  incident  and  the  scattering  grazing  angle,  respectively. 
(ps  is  the  histatic  angle.  Vww(8)  is  a  sound  reflection  coeffi¬ 
cient  from  the  seabottom  evaluated  at  the  grazing  angle  0.  ap 
is  the  complex  ratio  of  sediment  compressional  wave  speed 
to  water  sound  speed  (vp  is  a  real  parameter  of  the  speed 
ratio,  vp  =  cb/cw).  ap  is  the  sediment- water  density  ratio,  8p 
is  the  ratio  of  the  imaginary  and  real  parts  of  sound  speed  in 
the  sediment.  ocp  is  the  sound  attenuation  in  the  sediment  in 
dB/m. 


B.  The  SVS  model 


The  geometry  relevant  to  the  simplest  fluid  SVS  models 
is  shown  by  Fig.  14.1  of  Ref.  33.  If  the  volume  scattering 
strength  is  assumed  to  be  depth  independent  and  can  be 
treated  in  terms  of  an  effective  interface  scattering  strength 
(the  decibel  equivalent  of  scattering  cross  section  per  unit 
solid  angle  per  unit  area),  the  equivalent  interface  bistatic 
cross  section  for  SVS  can  be  expressed  by 


'Vv.s  {8j,  8S) 


\Vwp(di)\2\Vwp{ds)\2(TV 

2k>p|2Im[R(0i)+R(0s)]’ 


(22) 


<Jrbs(9i,9s)  =  kl\Aww\2W(U£), 

For  the  homogeneous  fluid  model, 

Awh-  =  ^  [1  +  K™,(0,')][1  +  VWW(9S)}G, 


(13) 

(14) 

(15) 


where 


P(9)  =  \Jctp2  —  cos2  0. 


(23) 


The  pressure  transmission  coefficient  ( Vwp )  is 

T/  _  WP 

wp~  z  +  1  ’ 


(24) 
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zwp(Q)  = 


apap  sin  8 , 
sin  6p 


(25) 


According  to  Jackson  and  Richardson,  Eqs.  (22)-(23)  “form 
the  basis  for  several  SVS  models,  which  differ  only  in  the 
assumptions  used  to  obtain  the  volume  scattering  cross  sec¬ 
tion,  oy.  The  simplest  approach  is  to  treat  av,  empirically, 
that  is  as  a  quantity  that  must  be  obtained  by  fitting  data 
rather  than  calculated  on  the  basis  of  theory.”1’  oy  has  the 
units  of  nr /m3  =  m  1 .  The  transmission  coefficient  Vwp  and 
the  seabottom  reflection  coefficient  Vww  have  a  relationship 

Vwp  =  1  +  K™,.  (26) 


Thus  Eq.  (22)  can  also  be  rewritten  as 


^svsiPh  Os ) 


\l  +  V„wm\2\l+V«U0s)\2vv 

2kw\ap\2lm[P(8i)  +  P(8s)} 


(27) 


IV.  INTEGRATING  THE  ENERGY  FLUX  METHOD  OF 
REVERBERATION  WITH  PHYSICS-BASED  SEABED 
SCATTERING  MODELS 

A.  In  the  angular  domain 

Equation  (12)  forms  a  general  SW  reverberation  expres¬ 
sion  in  the  angular  domain  with  the  physics-based  seabed 
scattering  model.  It  is  a  double-integral  over  incident/scatter¬ 
ing  grazing  angles.  All  classic  seabed  scattering  cross  sections 
expressed  as  a  function  of  incident/scattering  grazing  angles 
and  frequency  can  be  well  integrated  with  two-way  propaga¬ 
tion  in  the  angular  domain,  including  the  RBS  model  and  the 
SVS  model,  expressed  by  Eq.  (13)  or  (22). 

For  a  special  case,  in  the  Pekeris  waveguide  (iso-veloc- 
ity  water),  the  cycle  distance  of  the  nth  mode  ray  is  given  by 


pulse  duration,  to  (or  less  than  a  time  window  for  rever¬ 
beration  data  processing).  For  the  Pekeris  waveguide 
(iso-velocity  water),  Eq.  (29)  is  much  simpler  for  calcu¬ 
lating  SW  reverberation. 


B.  In  the  modal  domain 

From  Eqs.  (3)  and  (9),  we  have 
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(31) 


Using  Eqs.  (3),  (30),  and  (31),  the  SW  reverberation  expres¬ 
sion  that  integrates  over  mode-ray  angles,  Eq.  (12),  can  be 
converted  back  to  a  summation  over  normal  modes 


R(r,z-,zo)  =  ^-^-£E 


|0„(;o)|-  |Q„(z)r 

kmk„ 


x  e  2(A"+Wr  x  Acr(dm,8n) 
x  |0m(/t)|2|0„(/7)|2. 


(32) 


For  the  slowly  depth-varying  energy  dependence  of  the 
normal  mode  in  Eqs.  (3)  and  (32),  Ellis7  and  Yang6  had  the 
following  expression: 
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Sm  =  2/t/tan(0m)  -  l(dcph/dd)\g=gJkwsmdn} 
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Equation  (12)  becomes 
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To  overcome  a  problem  when  a  receiver  is  near  a  turning 
point,  where  kw(z )  «  kn(z),  alternatively,  Zhang  and  Jin 
have  proposed  the  following  expression:1 


l<WI2 


_ 2 K _ 

Sn[kl{z)D{z)  +  kl{z)~k^l2' 


(34) 


where 


In  derivation,  the  seabed  scattering  area  at  mid-  to  long- 
range  has  been  assumed  as  A  =  2 nr  x  dr  «  2 nr 
x(cH, to/2)  =  nrcwXQ.  The  waveguide  dispersion  spread¬ 
ing7  has  been  assumed  to  be  much  less  than  the  signal 


D(z)  =  0.875 


1  dc(z ) 
nf  dz 


2/3 
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Using  Eqs.  (8)  and  (34),  Eq.  (32)  can  be  re-written  as 
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(36) 
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This  is  a  simple  and  practical  expression  for  SW  reverbera¬ 
tion  in  the  modal  domain.  For  arbitrary  velocity  profiles,  Eq. 
(36)  can  be  used  to  calculate  SW  long-range  reverberation 
intensity  by  using  the  classic  cross  sections  of  seabottom 
scattering  and  any  available  normal-mode  code.  It  only 
requires  two  parameters:  Real  and  imaginary  parts  of  the 
modal  eigenvalue,  km  and  /)„, . 


are  both  in  the  water.  Substituting  Eq.  (41)  into  (40) 
results  in 


\^(h)\2 


i  +  vww{6»,)\2 ' 


(42) 


Thus  the  modal  scattering  matrix  in  SW  waveguides,  Eq. 
(38),  can  be  re-written  as 


V.  THE  MODAL  SCATTERING  MATRIX  DERIVED  FROM 
THE  ENERGY  FLUX  METHOD 

A.  A  simple  relationship  between  the  modal  scattering 
matrix  and  the  seabed  scattering  cross  sections 

Equation  (32)  can  be  re-written  as 


/?(r,z;z0)  =  ^e  53 


|O„,(z0)|2  |<F„(z)|2 


k  k 

x  e-2(ft»+A,)r  x  A  x  MSM(m,n). 


(37) 


Here  MSM(m ,  n)  is  defined  as  the  modal  scattering  matrix  in 
the  SW  waveguide  between  mode  m  and  mode  n,  expressed 
by  the  classic  seabed  scattering  cross  sections 

MSM(m,  n)  =  a(9m,  0„)  x  |4>m(/7)|2  |4>„(/2)|2.  (38) 


A  simple  relationship  between  the  waveguide  modal  scatter¬ 
ing  matrix  and  the  plane  wave  scattering  cross  section  is 
derived  as  follows. 

The  continuity  of  pressure  and  normal  component  parti¬ 
cle  velocity  at  the  water-sediment  interface  requires  that  the 
normal  mode  function  and  its  derivative  with  respect  to 
depth  z  have  the  relationship24'41 


qm  = 


1  d<S>„,  (z) 

0„,(z)  dz 


=  (Pw/Pb)y/%  ~  kl  (39) 


where  the  subscripts  w  and  b  denote  the  water  and  bottom  at 
the  interface  z  =  h.  Using  Eq.  (39),  Eq.  (33)  can  be  written  as 


MSM(m ,  n) 


|1  +  V„w(8m){ 


x  cr($w,  0n ) 


li  +  vw^f 


(43) 


Here  6m  =  cos-1  [km/kw (/?)].  $„,(/;)  is  the  amplitude  of  mth 
mode  at  the  water-bottom  interface.  o(()m.  ()„)  can  be  either 
the  RBS  cross  section  or  the  SVS  cross  section. 

Using  a  two-way  coupled-mode  representation,  Stotts 
and  Koch  first  derived  the  expression  of  Eq.  (43)  for  the 
RBS  model.29  Equation  (43),  derived  from  the  energy  flux 
method  for  SW  reverberation  in  this  paper,  shows  that  it 
should  be  applicable  to  any  seabed  scattering  cross  section, 
including  the  RBS  and  SVS  models. 

The  relationship  between  the  plane  wave  scattering 
cross  section  n(Qm.  9nj  and  the  waveguide  mode  scattering 
matrix  MSM,  expressed  by  Eq.  (43),  has  an  interesting  phys¬ 
ical  explanation.  9„)  is  the  classic  seabed  scattering 

cross  section  per  unit  seabed  area.  It  is  defined  as  the  scatter¬ 
ing  intensity  at  a  unit  distance  caused  by  an  incident  plane 
wave  with  unit  amplitude;  for  the  MSM  in  a  waveguide,  the 
equivalent  incident  wave  intensity  (decomposed  from  /??th 
mode)  and  the  scattering  intensity  (decomposed  from  nth 
mode)  at  a  unit  distance,  are  scaled  by  a  factor  of 
&m{h)/\l  +  Vmv(dm)\2  and  O2(/O/|l  +  Vw(0„)|2,  respec¬ 
tively.  The  factor  is  approximately  corresponded  to  the 
down-  or  up-going  (incident  or  scattering)  wave  intensity  at 
the  water-bottom  interface,  decomposed  from  a  mode  WKB 
expression. 
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Using  Eqs.  (18)  and  (24)-(26),  (or  using  8.46,  8.49,  and  8.53 
in  Ref.  33),  we  have 
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|1  +  Vww(6m)  |" 
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4 
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{  1  Pw 
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(41) 


Here  Vmv  is  the  reflection  coefficient  with  the  scripts 
ww  added  to  indicate  that  the  incident  and  reflected  fields 


B.  The  modal  scattering  matrixes  for  the  RBS  model 
and  SVS  model 

Substituting  Eqs.  (13)  and  ( 1 5)— (20)  into  Eq.  (43),  we 
have  the  MSM(m ,  «)for  the  RBSg  model  as  follows: 
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qn  is  defined  by  Eq.  (39).  The  MSM  for  the  RBS  model, 
derived  from  the  integration  of  the  reverberation  energy  flux 
method  with  RBS  model,  expressed  by  Eqs.  (44)  and  (45),  is 
same  as  the  result  derived  by  Tracey  and  Schmidt"  ’  and  by 
Song  et  al.,24  using  the  boundary  perturbation  method.  Song 
et  al.  called  bmn  as  “the  modal  backscattering  strength” 
between  modes  m  and  n.24 

In  the  same  way,  substituting  Eq.  (27)  into  Eq.  (43),  we 
have  the  MSM  for  the  SYS  model 


MSMSvs{m,n )  = 


1  +  Vmv(9m)\ 

IM0m)|2|M0B)|2gy 

K  2kw\ap\2lm[P(9m)  +  P(d„)} 
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(46) 


where  Vwp  is  the  pressure  transmission  coefficient  of  the 
compressional  wave  from  the  water  into  the  bottom.  p(9)  is 
defined  by  Eq.  (23). 


VI.  DATA-MODEL  COMPARISONS 

For  validating  the  energy  flux  method-derived  reverber¬ 
ation  expressions  integrated  with  the  HF  RBS  and  SVS  mod¬ 
els,  some  examples  are  shown  in  this  section. 


A.  Reverberation  measurements 

Reverberation  data  were  collected  on  June  3,  2001  at 
the  Asian  Sea  International  Experiment  (ASIAEX)  site  in 
the  East  China  Sea  using  38-  and  1000-g  explosive  charges. 
The  reverberation  levels  (RL)  averaged  over  three  shots 
were  analyzed  as  a  function  of  time  for  different  center  fre¬ 
quencies  in  1/3-octave  bands  and  normalized  to  the  source 
level  (SL).14  The  source  depths  were  around  50  m;  the  re¬ 
ceiver  depth  was  56.5  m.  The  horizontal  distances  between 
the  sources  and  the  receiver  were  in  a  range  of  8.2-15.7  m. 
Thus  the  RL  can  be  well  described  as  the  monostatic.  (More 
details  on  the  RL  measurements  can  be  found  in  Refs.  14 
and  44.)  The  water  depth  was  about  104  m.  The  sea  bottom 
was  nominally  flat.  During  the  reverberation  measurements 
on  June  3,  the  sea  surfaces  were  relatively  calm:  Wind  speed 
and  the  rms  surface-wave  height  were  2.74  m/s  and  0.10  m, 
respectively.44  The  average  sound  velocity  profile  at  the 
receiving  array  site  is  plotted  in  Fig.  1 . 


B.  Seabed  geoacoustic  and  scattering  models 

A  recent  paper  on  broadband  geoacoustic  inversions 
from  LF  field  measurements,  made  at  20  locations  around  the 
world  by  many  investigators,  has  resulted  in  an  effective 
seabed  geoacoustic  model.42  The  sound  attenuations  in  sandy 
or  sand-silt  mixture  bottoms,  inverted  from  different  acoustic 
field  characteristics,  exhibit  similar  magnitude  and  nonlinear 
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FIG.  1.  (Color  online)  Sound  velocity  profile  at  the  ASIAEX  site  during  the 
reverberation  measurements. 


frequency  dependence  at  all  of  these  locations,  including  the 
ASIAEX  site.  The  LF  field-derived  sound  speed  and  attenua¬ 
tion  in  the  sandy  bottoms  satisfy  the  Kramers-Kronig  relation¬ 
ship,4"  and  may  be  equally  well  described  by  four  physics- 
based  models:  The  Biot-Stoll  model,  the  Buckingham  VGS 
model,  the  Chotiros-Isakson  BICSQS  model,  and  the  Pierce- 
Carey  model.45  In  this  paper,  we  use  a  Biot  geoacoustic 
model  that  matches  the  LF  field  measurements.  All  the  effec¬ 
tive  Biot  parameters  are  listed  in  a  column  of  “LF  fit,”  Table 
10  of  Ref.  42.  Related  bottom  sound  attenuation  and  the  ve¬ 
locity  ratio  at  the  bottom-water  interface,  derived  from  the 
effective  Biot  model,  are  listed  in  Table  1. 

The  RBS  model  and  SVS  model  introduced  in  Sec.  Ill 
are  used  for  reverberation  computations. 


C.  Reverberation  data  are  in  good  agreement  with  the 
theoretical  model 

The  reverberation  expression  in  the  modal  domain, 
Eq.  (36),  is  used  to  calculate  RL  in  this  section.  First,  we 
assume  that  the  “spectral  exponent”  y2  of  the  seabed  roughness 
spectrum  in  Eq.  (14)  equals  3.0  (as  many  HF  seabed  scattering 
papers  suggested),  then  we  adjust  the  “spectral  strength”  CO2  to 
find  the  best  match  (by  vision)  between  the  numerical  RLs  and 
the  experimental  data.  For  the  SVS  model  in  Eq.  (27),  we 
adjust  the  parameter  oy  to  find  the  best  match.  The  data-model 
comparisons  at  800,  1500,  and  2500  Hz  are  shown  in  Figs. 
2(a)  to  Fig.  2(c).  These  figures  show  that  both  the  RBS  model 
and  SVS  model,  well  developed  by  the  HF  seabed  scattering 
community,  can  be  used  to  predict  LF  SW  reverberation  using 
a  set  of  right  seabed  geoacoustic  and  scattering  parameters. 
The  SW  reverberation  expression  derived  from  the  integration 
of  the  energy  flux  method  with  these  seabed  scattering  models 
is  also  validated  by  Figs.  2(a)-2(c),  as  well  as  by  the  results  in 
the  next  section. 

Figures  2(a)-2(c)  show  that  the  measured  reverberation 
levels  at  very  short  ranges  are  larger  than  the  theoretical  pre¬ 
dictions.  This  result  indicates  that  the  modal  expression  for 
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TABLE  I.  Bottom  attenuation  and  velocity  ratio  at  the  water-bottom  interface,  derived  from  the  Biot  model  that  matches  the  LF  field  measurements  (Ref.  42). 


/(Hz) 

150 

200 

300 

400 

500 

600 

700 

800 

1000 

1250 

1500 

1750 

2000 

2250 

2500 

c„/cw 

1.0604 

1.0605 

1.0606 

1.0608 

1.0611 

1.0614 

1.0618 

1.0622 

1.0632 

1.0646 

1.0663 

1.0681 

1 .0700 

1.0721 

1.0741 

dp  (dB/m) 

0.0110 

0.0187 

0.0401 

0.0694 

0.1064 

0.1509 

0.2026 

0.2613 

0.3979 

0.6012 

0.8344 

1.0909 

1.3643 

1.6486 

1.9386 

ASIAEX  site  1  =  800  Hz  SD  =  50  m  RD  =  56.5  m 


ASIAEX  site  f  =  1500  Hz  SD  =  50  m  RD  =  56.5  m 


ASIAEX  site  f=  2500  Hz  SD  =  S0m  RD  =  56.5  m 


FIG.  2.  (Color  online)  Comparisons  between  theoretical  RL(t)  (using  the 
RBS  and  SVS  models)  and  reverberation  data  at  given  frequencies,  (a)  800 
Hz,  (b)  1500  Hz,  (c)  2500  Hz. 
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SW  reverberation,  Eq.  (36),  can  not  be  used  for  predicting  the 
SW  reverberation  at  short  range.  One  of  the  reasons  is  that 
the  reverberation  modeling  in  this  paper  assumes  that  the 
amplitudes  of  the  down-  and  up-going  waves,  decomposed 
from  the  modes,  are  equal.  This  assumption  requires  that  the 
bottom  reflection  coefficient  is  close  to  one,  i.e.,  the  reverber¬ 
ation  data-model  comparison  should  be  made  at  long  ranges 
where  the  low-order  modes  with  low-grazing  angles  dominate 
the  sound  field.  In  fact,  the  modal  wave  numbers  are  complex 
(due  to  the  bottom  attenuation),  the  amplitudes  for  the  down¬ 
going/incident  wave  and  up-going/scattering  wave,  decom¬ 
posed  from  the  mode  WKB  expression,  are  not  equal.  This 
factor  can  not  be  ignored  for  short-range  reverberation. 

VII.  CHARICTERIZATION  OF  THE  LF  SEABOTTOM 
SCATTERING  FROM  LONG-RANGE  REVERBERATION 

A.  RL(t)  at  a  given  single  frequency  is  used  to  invert 
the  scattering  parameters  for  the  RBS  and  SVS  model 

Assuming  the  y2  =  3.0,  the  reverberation  level  as  a 
function  of  time,  RL(t),  at  given  frequencies  is  used  to  find 
the  best  match  between  measurements  and  models,  as  we  did 
in  Sec.  VI.  The  theoretical  RL  is  numerically  evaluated  by 
using  Eq.  (36)  and  a  mode  code  of  Modal  Acoustic 
Transmission  Loss  (MOATL),46  using  different  cot  for  the 
RBS  model  or  ay  for  the  SVS  model.  When  the  numerical 
RL  curve  matches  the  experimental  RL(t)  data,  co 2  or  oy  is 
determined  at  the  given  frequency.  All  resultant  scattering 
parameters  from  the  best  matches  for  both  the  RBS  and  SVS 
models  at  different  frequencies  are  listed  in  Table  II. 

1.  Internal  parameter  coupling  between  y2  and  o>2  in 
the  RBS  model 

For  the  HF,  the  bottom  roughness  is  assumed  to  be  a 
random  process  described  by  (a  two  parameter,  isotropic, 
two  dimensional  roughness  spectral  density)  Eq.  (14).  In 
general,  y2  and  co2  are  assumed  as  constant.  The  roughness 
spectrum  exponent  is  restricted  to  the  range  of  2  <  72  <  4 
with  a  mean  of  3.0.  Figures  2(a)-2(c)  and  Table  II  show  that 
when  the  spectral  exponent  y2  =  3.0,  the  data-model  fittings 
at  different  frequencies  require  the  RBS  model  to  have  very 
different  values  of  the  spectral  strength  m2.  Numerical  simu¬ 
lations  show  that  at  a  given  single  frequency,  there  are  vari¬ 
ous  combinations  of  y2  and  co2  in  the  RBS  mode  that  can 
match  the  measurements.  As  an  example,  Figs.  3(a)  and  3(b) 
show  that  three  different  pairs  of  y2  and  co2  would  result  in 
an  identical  RL(t).  An  increase  (decrease)  in  y2  can  be  com¬ 
pensated  for  by  an  increase  (decrease)  in  012-  Thus  it  is 
impossible  to  uniquely  invert  the  bottom  roughness  spectral 
density  from  SW  reverberation  at  one  single  frequency.  It 
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TABLE  II.  RL(t) -inverted  LF  seabed  scattering  parameters,  assuming  y2  —  3.0. 


/(Hz) 

150 

200 

300 

400 

500 

600 

700 

800 

1000 

1250 

1500 

1750 

2000 

2250 

2500 

C02  (  X  1 0  4) 

0.6 

0.8 

0.9 

1.2 

2.7 

3.1 

4.6 

5.4 

8.0 

18.0 

24.0 

28.0 

34.0 

42.0 

34.0 

cy(x  10~4) 

0.029 

0.067 

0.16 

0.43 

1.50 

2.42 

5.68 

8.90 

20.0 

72.0 

142.0 

240.0 

382.0 

626.0 

620.0 

needs  wideband  reverberation  measurements  as  described  in 
Sec.  VII B. 


2.  Strong  frequency  dependence  of  the  LF  volume 
scattering  cross  section  av 

In  the  HF  range,  the  sediment  volume  cross  section, 
ay,  is  often  quantified  by  the  following  dimensionless 
parameter33 

(72  =  av/aip,  or  av  =  a2Xp,  (47) 

where  y.p  is  the  bottom  attenuation  in  dB/m.  For  high  fre¬ 
quencies,  ap  is  close  to  linear  frequency  dependence,  a2  is 
taken  to  be  frequency  independent,  then  Eq.  (22)  or  (27)  will 
yield  a  scattering  cross  section  (erm)  that  is  nearly  frequency 


ASIAEX  site  f  =  1500  Hz  SD  =  50/77  RD  -  56.5  m 


ASIAEX  site  f=  2500  Hz  SD  =  SO  m  RD  =  56.5  m 


FIG.  3.  (Color  online)  Three  different  pairs  of  y2  and  oj2  result  in  an  identi¬ 
cal  RL(t).  An  increase  (decrease)  in  y2  can  be  compensated  for  by  an 
increase  (decrease)  in  co2. 

62  J.  Acoust.  Soc.  Am.,  Vol.  1 34,  No.  1 ,  July  201 3 


independent.33  However,  Figs.  2(a)-2(c)  and  Table  II  show 
that  at  low  frequencies,  the  reverberation-derived  volume 
scattering  cross  section,  ay,  has  much  stronger  frequency 
dependence. 

All  ay  in  Table  II,  inverted  from  RL(t)  at  different 
frequencies,  are  plotted  in  Figs.  4(a)  and  4(b)  against  the 
bottom  attenuation  and  frequency,  respectively.  Data  fitting 
results  in  the  following  approximate  relations: 

av  «  O.O17O0C3057,  (48) 

av  =  (0.00235  ±0.000 19)(f/l  000)  (3'822±0092) ;  (49) 

where  a.p  is  the  bottom  attenuation  in  dB/m, /is  frequency  in 
Hz. 

B.  Wideband  RL(f)  at  giving  reverberation  times  is 
used  to  invert  the  scattering  parameters  for  the  RBS 
and  SVS  models 

The  reverberation  level  as  a  function  of  frequency, 
RL(f),  at  given  times  is  used  to  invert  the  LF  seabottom  scat¬ 
tering  at  low  grazing  angles.  For  the  RBS  model,  both  y2  and 
co2  are  assumed  as  adjustable.  For  the  SVS  model,  as  men¬ 
tioned  in  the  last  section,  the  HF  expression  of  Eq.  (47) 
needs  to  be  modified  for  LF  reverberation  modeling.  We 
assume 

av  =  0*2*7  ■  (50) 

For  a  given  time,  a  theoretical  RL  as  a  function  of  frequency, 
RL(f),  is  numerically  evaluated  by  using  Eq.  (36)  and  a 
mode  code  of  MOATL  for  a  pair  of  y2  and  a>2  (for  the  RBS 
model)  or  a  pair  of  a2  and  riy  (for  the  SVS  model).  The  pair 
of  y2  and  c o2,  or  the  pair  of  a*  and  ny,  for  which  the  numeri¬ 
cal  RL(f)  curves  best  match  the  experimental  RL(f)  data,  is 
both  readily  and  uniquely  determined.  The  resultant  best 
matches  between  theoretical  predictions  and  measurements 
are  shown  in  Figs.  5(a)  to  5(h)  All  the  AT  (/'/inverted  param¬ 
eters  at  different  times  for  the  bottom  roughness  spectrum 
(y2  and  <n2)  and  for  the  SVS  cross  section  [a*  and  ny  in 
Eq.  (50)]  are  listed  in  Table  III.  The  standard  deviations 
between  the  measurements  and  theoretical  results  are  also 
listed  in  Table  III. 

C.  Characteristics  of  the  LF  RBS  model  and  SVS 
model  at  low  grazing  angles 

To  ensure  that  the  characteristics  of  seabed  scattering 
derived  from  reverberation  are  limited  to  low  frequencies 
and  low  grazing  angles  and  to  reduce  possible  bottom  scat¬ 
tering  contributions  from  deeper  layers  of  sediments,  in  this 
section,  we  just  use  the  results  derived  from  RL(f)  in  a  range 
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from  RL(t)  at  different  frequencies 


Frequency  (  Hz  ) 


FIG.  4.  (Color  online)  The  sediment  volume  scattering  cross  section  ay, 
inverted  from  RL(t)  at  different  frequencies,  (a)  Gy  vs  the  bottom  attenua¬ 
tion,  (b)  Gy  vs  frequency. 


of  3-12  s;  the  ranges  from  the  source  to  the  seabed  scattering 
areas  are  about  22-88  times  of  the  water  depth.  The  choice 
of  the  RL  data  between  3  and  12  s  is  based  on  the  following 
considerations:  (1)  As  we  discussed  in  a  previous  paper,14  to 
get  quality  wideband  reverberation  levels  (RL)  normalized 
to  the  source  level  (SL)  simultaneously  for  both  short  and 
long  range  is  a  very  delicate  task  that  can  be  subject  to 
errors.  At  short-range,  the  measurements  often  encounter 
signal  overflow  or  saturation  and  possible  nonlinear  gain 
problems  in  the  reverberation  measurement  system  caused 
by  a  powerful  explosive  source.  (2)  Due  to  strong  environ¬ 
mental  noise  interference,  it  is  very  hard  to  get  long-range 
RL  with  a  higher  reverberation-to-noise  for  lower  frequen¬ 
cies,  such  as  150  and  200  Hz.  Thus  the  quality  of  RL  data 
between  3  and  12  s  should  be  more  reliable  when  a  broad¬ 
band  of  150-2500  Hz  is  used.  (3)  The  RBS  cross  section  arbs 
and  the  equivalent  interface  bistatic  cross  section  for  SVS 
asvs,  expressed  by  Eqs.  (13)  and  (27),  are  derived  from  the 
perturbation  approximations.  They  are  more  accurate  at  low- 
to  mid-grazing  angles,  and  not  good  for  very  large  grazing 
angles.  The  RL  data  between  3  and  12  s  ensure  all  inverted 
seabed  scattering  are  for  low  grazing  angles  in  the  whole  fre¬ 
quency  band  of  150-2500  Hz.  (4)  The  reverberation 


expressions  in  this  paper  are  better  for  long  range  reverbera¬ 
tion  predictions  (see  the  discussion  made  in  the  last  para¬ 
graph  of  Sec.  VI).  (5)  The  seabed  effective  geoacoustic 
model  in  this  paper  is  derived  from  LF  long-range  held 
measurements  and  is  assumed  as  a  homogenous  half  space. 
The  long-range  reverberation  data  ensure  only  a  top  layer  of 
the  bottom  is  under  the  consideration  for  the  LF  seabed  scat¬ 
tering  characterization. 

1.  The  LF  bottom  roughness  spectrum  from 
long-range  reverberation 

In  contrast  to  Figs.  2  and  3,  obtained  from  RL(t)  at  given 
single  frequencies,  Fig.  5  and  Table  III  show  that  the  wide¬ 
band  reverberation  data  as  a  function  of  frequency,  RL(f), 
between  150-2500  Hz  at  different  reverberation  times  may 
uniquely  determine  the  bottom  roughness  spectrum.  For 
example,  based  on  Table  III,  the  average  LF  bottom  rough¬ 
ness  spectrum  parameters  for  low  grazing  angles  at  the 
ASIAEX  site,  derived  from  RL(f)  data  between  3  and  12  s  in 
a  frequency  band  of  150-2500  Hz,  are  close  to 

y2  ~  1-3,  co2  «  2.8  x  10~5.  (51) 

These  values  are  different  from  those  commonly  used  for  HF. 

Table  III  also  shows  that  with  decreasing  time/range 
(corresponding  to  increased  grazing  angles  of  the  seabed 
scattering),  RL(f) -inverted  y2  and  co2  seem  to  increase.  y2, 
inverted  from  short-range  reverberation  at  1.5  and  2  s,  may 
reach  to  2. 2-2. 3.  The  reason  for  this  is  not  clear  at  this 
moment.  It  might  be  due  to  the  scattering  contributions  from 
deeper  layers  of  sediments  at  short  range  for  lower  frequen¬ 
cies  or  due  to  the  problems/physics  discussed  in  the  first 
paragraph  of  Sec.  VII C.  In  addition,  the  theoretical  RL  val¬ 
ues  at  shorter  range  are  also  less  accurate.  Thus  the  results 
from  1.5  and  2.0  s  are  included  in  Table  III  just  for  refer¬ 
ence/discussion. 

2.  The  LF  volume  scattering  cross  section  nv  from 
long-range  reverberation 

All  the  RL(f)  inverted  cy  from  3  to  12  s  in  Table  III  are 
plotted  in  Fig.  6(a).  Using  the  power  law  fitting,  the  average 
volume  cross  section  in  sediments  can  be  expressed  by 

crv  =  0.0171a),2  038±0001).  (52) 

Using  the  frequency  dependence  of  the  bottom  attenuation 
listed  in  Table  I,  all  the  ay  obtained  from  3  to  12  s  are  plotted 
against  the  frequency  in  Fig.  6(b).  From  the  power  law  fit¬ 
ting,  the  average  volume  cross  section  from  all  RL(f) 
between  3  and  12  s  can  be  expressed  by 

gv  =  (0.00238  ±0.00004)  (// 100)  (3'805±0021).  (53) 

Equations  (52)  and  (53)  are  very  close  to  Eqs.  (48)  and  (49), 
where  txp  is  the  bottom  attenuation  in  dB/m,/is  frequency  in 
Hz.  The  long-range  reverberation-derived  LF  oy  exhibits 
strong  frequency  dependence  that  is  very  different  from  that 
in  the  HF  range.  The  results  in  Fig.  6  are  also  almost  the 
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ASIAEX  Site  150  -2500  Hz.  t  =  3.0  (sec),  r/h~22 


ASIAEX  Site  150 -2500  Hz,  t=  7.0  (sec),  r/h~51 


ASIAEX  Site  ISO-  2500  Hz,  t  =  4.0  (sec),  r/h~29 
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FIG.  5.  Comparisons  between  theo¬ 
retical  RL(f)  (using  the  RBS  and 
SVS  models)  and  reverberation  data 
at  given  times,  (a)  3.0  s,  (b)  4.0  s,  (c) 
5.0  s,  (d)  6.0  s,  (e)  7.0  s,  (f)  8.0  s,  (g) 
10.0  s,  and  (h)  12.0  s. 


same  as  those  in  Fig.  4,  obtained  from  RL(t)  at  different 
frequencies. 

As  this  time,  we  do  not  know  why  the  LF  bottom  rough¬ 
ness  spectrum  and  the  SVS  cross  section  at  low  grazing 
angles,  derived  from  wideband  long-range  reverberation, 
have  different  characteristics  from  those  commonly  used  for 
HF.  We  also  do  not  know  which  of  these  should  be  a 


dominant  LF  seabed  scattering  mechanism  for  long-range 
reverberation.  These  questions  are  outside  of  the  present 
paper’s  scope.  Hopefully,  the  LF  seabed  scattering  parame¬ 
ters  at  low  grazing  angles  inverted  from  long-range  reverber¬ 
ation  measurements  in  this  paper  can  offer  some  reference 
data  sets  for  future  analysis  of  the  LF  seabed  scattering 
mechanisms. 
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TABLE  III.  LF  seabed  scattering  parameters  inverted  from  RL(f)  at  different  reverberation  times. 


Models 

RBS: 

Bottom  roughness  spectrum  W(AK) 

=  m/(\K)y2 

SVS:  Volume  scattering  cross  section:  oy  — 

Rev.  time  (s) 

r/hn 

12 

0)2  X  10-5 

StdRBs 

°2 

nv 

Stdsvs  (dB) 

1.5* 

11 

2.3 

23.7 

2.031 

0.010 

1.5 

2.294 

2.0* 

15 

2.2 

18.3 

1.411 

0.009 

1.5 

1.649 

3.0 

22 

1.6 

4.9 

1.419 

0.014 

1.9 

1.603 

4.0 

29 

1.4 

3.2 

1.192 

0.015 

2.0 

1.348 

5.0 

37 

1.3 

2.6 

1.719 
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VIII.  SUMMATION  AND  DISCUSSION 

The  energy  flux  (angular  spectrum)  method  for  SW 
reverberation,  based  on  the  WKB  approximation  to  the 
normal-mode  solution,  has  been  integrated  with  the  physics- 
based  RBS  and  SVS  models.  The  integration  also  results  in  a 
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FIG.  6.  (Color  online)  The  sediment  volume  scattering  cross  sectionay, 
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simple  relationship  between  the  classic  seabed  scattering 
cross  sections  and  the  modal  scattering  matrix  in  wave¬ 
guides.  Taking  advantage  of  the  progress  made  by  the  seabed 
scattering  community,  and  using  commonly  known  seabed 
geoacoustic/scattering  models,  the  resultant  SW  reverbera¬ 
tion  expressions  in  the  angular  domain  or  in  the  modal 
domain  can  conveniently  be  used  for  predicting  SW  rever¬ 
beration  for  the  Pekeris  waveguide  or  non-Pekeris  wave¬ 
guide,  respectively.  Comparisons  of  the  reverberation  model 
with  the  reverberation  data,  collected  from  the  ASIAEX  site 
in  a  frequency  range  of  150-2500  Hz,  show  that  the  the  RBS 
model  or  the  SVS  model  can  directly  be  used  for  the  LF 
reverberation  modeling. 

The  reverberation  data  at  the  ASIAEX  site  are  used  for 
characterization  of  the  LF  seabed  scattering  at  low  grazing 
angles.  The  results  show  that  the  reverberation  level  as  a 
function  of  time,  RL(t),  at  one  single  frequency  cannot 
uniquely  determine  the  seabed  scattering  parameters.  The 
wideband  reverberation  data  can  uniquely  determine  a  set 
of  the  bottom  roughness  spectra  or  the  SVS  cross  section. 
The  bottom  roughness  spectrum  and  the  SVS  cross  section, 
inverted  from  the  broadband  long-range  reverberation  data 
at  the  ASIAEX  site,  exhibit  different  characteristics  from 
those  commonly  used  in  HF.  For  example,  for  the  HF,  the 
roughness  spectrum  exponent  is  restricted  to  the  range  of 
2  <  y2  <  4  with  a  mean  value  of  3.0.  But  the  long-range 
reverberation  inverted  LF  roughness  spectrum  exponent 
y2  ~  1.3.  The  HF  SVS  cross  section  ay  is  generally 
assumed  to  have  linear  frequency  dependence.  However,  in 
the  150-2500  Hz  range,  the  long-range  reverberation- 
inverted  ay  exhibits  much  stronger  frequency  dependence: 
av  oc/3-8. 

All  the  reverberation  inversions  in  this  paper  are  con¬ 
strained  by  the  Biot  geoacoustic  model  for  the  bottom.  The 
reverberation  data  from  the  ASIAEX  site  can  almost  be 
equally  well  predicted  by  the  RBS  model  and  SVS  model. 
To  analyze  which  one  (or  both)  is  a  dominating  LF  scatter¬ 
ing  mechanism  at  low  grazing  angles,  it  is  desirable  to  have 
geophysical  measurements  on  the  bottom  roughness  spec¬ 
trum  and  sediment  inhomogeneities  for  the  seabed  scatter¬ 
ing  modeling.  Further  research  is  called  for  understanding 
the  mechanism  of  the  reverberation-inverted  LF  seabed 
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scattering  characteristics  and  analyzing  the  uncertainties  in 
reverberation  inversion  with  more  measurements. 
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